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Abstract 



The work in this thesis consists of two distinct parts: 

A class of models called, "String-flip potential models," (SFP's) are studied as 
a possible candidate for modeling nuclear matter in terms of constituent quarks. 
These models are inspired from lattice quantum-chromodynamics (QCD) and are 
nonperturbative in nature. It is shown that they are viable candidates for modeling 
nuclear matter since they reproduce most of the bulk properties except for nuclear 
binding. Their properties are studied in nuclear and mesonic matter. A new class 
of models is developed, called "flux-bubble potential models," which allows for the 
SFP's to be extended to include perturbative QCD interactions. Attempts to obtain 
nuclear binding is not successful, but valuable insight was gained towards possible 
future directions to pursue. 

The possibility of studying Superstring inspired E 6 phenomenology at high energy 
hadron colliders is investigated. The production of heavy lepton pairs via a gluon- 
gluon fusion mechanism is discussed. An enhancement in the parton level cross- 
section is expected due to the heavy (s)fermion loops which couple to the gluons. 
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Chapter 1 
Introduction 



Throughout the history of the universe many forces have played a role in shaping 
it into what it is today 0, ^j. In this thesis the role of Quantum Chromodynamics 
(QCD) in the creation of baryonic matter and the possible low energy consequences 
of the exotic theory of superstrings will be investigated. 

In particular, the viability of a class of models for nuclear matter called string- 
flip potential models will be considered. String-flip potential models attempt to 
explain nuclear matter from the more fundamental constituent-quark-level picture. 
Explaining nuclear matter in this way is by no means an easy feat: attempts to do 
so have met with varied results. Unlike nuclear physics, which is quite successful at 
explaining nuclear phenomena from a nucleon perspective, constituent quark models 
of nuclear matter are far from complete. The main stumbling block is the non- 
perturbative and many-body nature of the strong interaction. Very recently some 
inroads have been made in the area of lattice QCD that may soon prove to be 
revolutionary to this field |§. In fact the models that will be examined here were 
inspired by lattice QCD. 

The phenomenology of superstring-inspired E 6 models will also be investigated. 



1 



CHAPTER 1. INTRODUCTION 



2 



In particular, the possibility of heavy lepton production at high energy hadron collid- 
ers will be studied. Such a find would help solve the generational hierarchy problem 
in the standard model and lend support to a theory that unifies all the known forces 
of nature, namely superstrings. 



1.1 The String-Flip Potential Model 

For the past 30 years several attempts have been made, with very little success, to 
describe nuclear matter in terms of its constituent quarks. The main difficulty is 
due to the non-perturbative nature of QCD. The most rigorous method for handling 
multiquark systems to date is lattice QCD. However, lattice QCD is very compu- 
tationally intensive and given the magnitude of the problem it appears unlikely to 
be useful in the near future.0 As a result, more phenomenological means must be 
considered. 

The idea of string-flip potential models is borrowed from certain results in lattice 
QCD and experimental particle physics. A potential derived from computations in 
lattice QCD is confirmed by fitting mesonic spectra in particle physics experiments. 
It has been found that the most consistent inter-quark potential model between 
quarks and antiquarks has the form 

4 a 

V(r) ~ or --— . (1.1) 

3 r 

This formula is basically an interpolation between the long range non-perturbative 
(err) and short range perturbative (— | ^f-) parts of the force between pairs of quarks 



(figure |1 . 1| ) . The string-flip potential model ignores the short range part of the 



1 Some very recent advancements have been made in the area of lattice QCD that have reduced 
computation time by several orders of magnitude. "Now what took hundreds of Cray Supercom- 
puter hours can be done in only a few hours on a laptop computer." pfl. 
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potential and considers an ensemble of quark-antiquark pairs, qq , such that the 
total amount of string, J2 r qq > shared between them is minimized. 




(a) (b) 

Figure 1.1: The colour field lines between quarks collapse upon themselves, due to the self 
interacting nature of the gluons, to form a flux-tube-like structure. At long distances (a) the 
fields lines collapse to become almost string-like and at short distances (b) the fields lines 
expand to become almost QED-like. 

This particular model has been used in an attempt to model nuclear matter. 
Although it has an obvious shortcoming, in that it is more applicable to a pion gas, 
it does surprisingly well at predicting some of the overall bulk properties of nuclear 
matter (§ §. 

It is fairly straightforward to generalize this simple model to a more realistic 
one which involves triplets of quarks. Here the flux-tubes leaving each quark meets 
at a central vertex such that overall length, r, of "flux-tubing" is minimized (cf. 
figure |2.1| .a). The potential energy is simply or ||. 

In a more general setting one could consider a full many-body quark potential in 
which large clusters of quarks may be connected by a single network of flux-tubes ||; 
in general, this "gas" is assumed to consist of colourless objects. Again the potential 
energy is simply ar , where r is now the minimal amount of flux-tubing used for a 
given cluster of quarks. In such a model, there could exist very complex topological 
configurations of flux-tubes, such as long strands or web-like structures (figure |1.2| ). 
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Figure 1.2: Possible flux-tube configurations for an SU C (3) quark gas in a box. 

All of these models are completely motivated by results from lattice QCD, where 
variations are taken about minimal lattice field configurations between quarks. 



1.1.1 Possible Phases of Nuclear Matter 

There is a major advantage to understanding nuclear matter in terms of its con- 
stituent quarks. Not only is a deeper understanding of nuclear physics likely to be 
achieved, but also a more general understanding of the nature of quark matter. This 
understanding could possibly lead to the prediction of more exotic forms of mat- 
ter. To illustrate this point let us now "hypothesize" some of the possible phases 
of nuclear/quark matter (figure |1.3|) in the context of a general string-flip potential 
picture. 

First consider nuclear matter in a box at low temperature (i.e., in its ground 
state) in a standard nuclear model, with no references to quarks. At very low density 
the system is a gas of isolated nucleons with a Fermi degeneracy pressure on the walls 
of the box. As the box is gradually squeezed some of the nucleons may come close 
enough to start clumping together. At this stage the pressure becomes negative due 
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Figure 1.3: Postulated nuclear matter phase diagram {adapted from Siemens & Jensen /Q/. 
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to the clumping forces which are trying to reduce the volume. As the box is squeezed 
even further, short-range repulsion and other saturation mechanisms Jl|, |7j cause the 
pressure to become positive. Effectively the system behaves very much like water 
vapor; little droplets of nucleons floating around in a box. Further squeezing of the 
box causes liquid (probably a superfluid) to be formed. 

Now consider a very simple flux-tube model where the many-body potential 
contains only three-body forces between red (r), green (g), and blue (b) quarks: 
each nucleon is represented by three valence quarks, rgb, connected by a triangular 
web of flux-tubes (cf. figure |2.1|.a). In a flux-tube picture it is suspected that the role 
of meson exchange between nucleons in the nuclear physics picture is mimicked by 
the swapping of flux-tubes as nucleons move close to each other M. At sufficiently 
low pressure the system behaves as a gas of nucleons because the clusters of quarks 
are essentially isolated from one another: no flux-tubes are exchanged. In addition, 
as the quarks are fermions a Fermi pressure is set up on the walls of the box but 
not within the rgb clusters themselves. As the box is squeezed some of the clusters 
come into contact with each other causing a clumping effect. The saturation effect 
is perhaps not that obvious, however it has been shown that this simple flux-tube 
picture does indeed lead to saturation of nuclear forces [§], ||, ||. It even produces 
the subtle effect of nucleon swelling in nuclear matter.! However it does not appear 
to produce a strong enough attractive force to produce binding; speculation as to 
"why?" will be discussed later on in this thesis. "Further squeezing of the box causes 
a liquid (probably a superfluid) to be formed." Presumably if spin correlations were 
set up between pairs of clusters of quarks, collective states that indicate superfluidity 
may be detected. This idea has not been tested. The possibility of new physics is 
introduced when the box is squeezed further. At such densities the nucleon could 

2 The EMC effect 
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perhaps form a liquid crystal flQfl , as it might be more energetically favourable for 
the planes of the rgb clusters to align themselves. If the box is squeezed even further 
the flux-tubes essentially dissolve leaving a Fermi gas of quarks. 

If the temperature of the box is now increased, even more phases of nuclear 
matter become evident. For very low densities the system remains a Fermi gas of 
nucleons, but now it may become an excited Fermi gas. As the box is squeezed the 
nucleons may clump to form a gas-liquid mixture provided the temperature is not 
too high. If the temperature is too high it may simply stay in its gaseous phase. As 
the box is squeezed further each of the aforementioned possible states will go over 
into a liquid. If the temperature and pressure is just right it is possible for all three 
phases of nuclear matter to coexist: cf. figure P~7B| . 

If the temperature at low density is pushed higher the gas pionizes (figure |1.3|) : 
the mesons go on shell. In order for this to work in a flux-tube model a meson 
production mechanism would have to be incorporated. Perhaps the simplest incor- 



poration would be to introduce a string breaking mechanism, figure |1.4| , then if the 
flux-tubes became too long they would break, producing mesons. If the temperature 
is further increased a plasma of quarks, antiquarks, and gluons would be produced. 

Finally, by varying the temperature and pressure of the aforementioned phases 
other phases can be reached: if the temperature is increased in the liquid crystal 
phase (see figure |1.3p the liquid crystal will eventually become disrupted forming, 
perhaps, a pionized gas or a plasma. If the pressure or temperature of the liquid 
phase is raised a liquid crystal, pionized gas, or plasma may be formed (this suggests 
perhaps another triple point on the phase diagram). 

A summary of the various phases of nuclear matter that have been discussed 
is illustrated in figure It should be emphasized that all of this is specula- 

tive and model dependent, i.e., the only experiments that exist are in the region 
T < Q(l)MeV and p < 0(0.2) fni' 3 f§. 
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Figure 1.4: A naive meson production scheme. If a flux-tube is stretched long enough such 
that there is enough mass energy in a segment to produce a qq pair of mass 2m q then the string 
(flux-tube) will break creating a pair of quarks, with the corresponding segment vanishing, with 
probability ~ e ^ 7!m q / 17 pfl . Note that this model does not include a mechanism for getting rid 
of the mesons. 

1.1.2 A Crude Model of Nuclear Matter 

It should now be apparent that attempting to construct a model of nuclear matter 
in terms of quarks is not a very simple task. To reduce this burden some simplifying 
assumptions or restrictions will have to be made. One restriction is to only consider 
low temperature phenomena. Some simplifying assumptions would be to require 
that the model correctly predict some of the very basic bulk properties of nuclear 
matter: 

• nucleon gas at low densities with no van der Waals forces 

• nucleon binding at higher densities 

• nucleon swelling and saturation of nuclear forces with increasing density 

• quark gas at extremely high densities 

There are many models that attempt to reproduce all of these properties but none 
that reproduces them completely. 
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In this thesis only one particular class of models, called string-flip potential mod- 
els m [|, U, will be considered, in chapters 2 and 3. These models appear to be 
promising because they reproduce most of the aforementioned properties with the 
exception of nucleon binding. 



1.2 Superstring Inspired Eq Models 

The SU(3) C ® SU(2) L <g> U(l)y standard model (SM) is a very successful model. It 
has thus far withstood rigorous experimental testing. However, despite its success 
the SM has many problems: 

• no unification of the forces 

• gauge hierarchical and fine tuning problems 

• three generations of quarks and leptons for no particular reason 

• too many parameters to be extracted from experiment 

Some of the earlier attempts at unification tried to unify the strong and elec- 
troweak forces by embedding the SU(3) C £g> SU(2) L ® U(l)y structure into higher 
groups, such as SU(5) and SO(10). These "grand unified theories", or GUT's, were 
only partially successful. The simplest of the GUT's was SU(5) which seemed promis- 
ing at the time because it predicted the ratio of the SU(2) W and U(l) em couplings 



and the proton lifetime [12|. However, the ordinary SU(5) GUT is no longer a pos- 
sibility because more refined experimental measurements are now in disagreement 
with its predictions for the couplings and the proton lifetime |13[]. In addition, this 
simple model had too many parameters and no explanation for family replication. 
The next likely candidate group was SO(10), although the three (or more) copies of 
the generational structure still had to be inserted by hand. 
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Difficulties with the SM and GUT models concerning gauge hierarchy and fine 
tuning problems led to theoretical remedies such as technicolour and supersymmetry 
(SUSY). The most appealing of these theories is SUSY, which has generators that 
relate particles of different spin in the same sup er mult iplet. The locality of these 
generators leads to supergravity models. SUSY (and its extended versions) however, 
did not have enough room for all of the SM particles |12|]. To solve this problem 
direct product structures were made with SUSY and Yang-Mills gauge groups. These 
structures are now commonly referred to as "SUSY" models. Of course the price 
paid for this was a large particle spectrum (at least twice that of the SM) and the 
problem of family replication still remained. 

In the early 1970's some interest was sparked in E 6 as a GUT when it was 
discovered that all the then known generations of fermions could be placed in a single 
27 dimensional representation. This ("topless") model was quite popular because 
the newly discovered r lepton and b quark could also be fitted neatly into the 27; 
there was no need for a third generation. However this model was quickly disallowed 
as it was experimentally shown that the r and b belonged to a third generation, and 
the idea of E 6 as a GUT died. 

In late 1984 Green and Schwarz showed that 10 dimensional string theory is 
anomaly free if its gauge group is either Eg <g> E' 8 or SO (32). The group that had 
received the most attention was E 8 (g) E' 8 as it led to chiral fermions, similar to those 
in the SM, whereas SO(32) did not. Furthermore, it was shown that compactifica- 
tion down to 4 dimensions (assuming N=l SUSY) can lead to Eq as an "effective" 
GUT group. Each family of SM particles now sits in its own 27, figure |1.5| . The 
generational problem may be solved because it is expected that any reasonable com- 
pactification scheme should generate the appropriate number of copies of the 27. 
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Figure 1.5: Eg particle content. The SM particles are shown in the boxes on the left and 
their "exotic" counterparts outside the boxes on the right. Although the exotics are labeled 
in a way that suggests they have the same quantum numbers as the non-exotics, in general 
they need not. The labeling for these particles in the literature has not been settled upon and 
varies quite significantly from paper to paper Jl3[ . Here the labeling scheme was chosen to 
reflect a specific Eg model that will be constructed in this thesis. In particular, all the exotics 
will carry the "expected" quantum numbers as their non-exotic counter parts do, with the 
exception being L=0 (i.e., Lcpton number zero) for the primed and double primed ones. 



"Note: Embedded in the 27's is the symmetry group SU(2)i due to an ambiguity in the 
particle assignments \ ( v ] ) d£ ^ ^ \ ( "i ) d'° L and {if^} {v['^} (cf. figure O.d). This 



ambiguity can easily be seen via the decomposition 27=^ e (SO(10),SU(5)). 
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For instance in a Calabi-Yau compactification scheme [[14 

E 8 <g> Eg — ► SU(3) <g> E 6 <g> E' 8 , 

the number of generations is related to the topology of the compactified space. 
A further assertion that the matter fields remain supersymmetrically degenerate 
ensures proper management of any gauge hierarchical and fine tuning problems. It 
is assumed that the hidden sector, Eg, which couples to the matter fields of E 6 by 
gravitational interactions will provide a mechanism for lifting the degeneracy. 

So the inspiration for using E 6 is that if it proves to be a possible GUT then 
it opens up the possibility of finding a TOE (Theory Of Everything). However, it 
should be pointed out that E 6 is not the only possible stop en route to the SM, but 



it is the most studied |L3|. It is for this reason that the low energy phenomenology 



resulting from E 6 will be studied in this thesis. 

1.2.1 Eg Phenomenology 

An extra Ze 

In order to produce SM phenomenology E 6 must be broken. Also to handle any 
hierarchical and fine tuning problems, SUSY must be preserved Jl4[]. This restriction 



makes the task more difficult, using most naive breaking schemes. The solution to 



the problem was found by using a Wilson-loop mechanism |14| over the non-simply 



connected-compactified-string-manifold to factor out the various subgroups of E 6 . 



Figure [L6] shows some of the possible, popular, rank 5 and rank 6 groups that can 
be produced by this scheme. As it can be seen, the various breaking schemes 
always give rise to extra vector bosons beyond the SM: in fact it is unavoidable. In 
this thesis only the simplest of these models (figure |1.6| .a) which generates an extra 
vector boson, the Ze , will be considered. 



CHAPTER 1. INTRODUCTION 



13 



(a) E 6 — > SU(3) C ® SU(2) L <g> U(1) Y ®U(1)y b 

(b) E 6 ^SO(10)®U(1)^ ) 

I } SU(3) C ® SU(2) L ® U(1) Y ®U(l) fl 

SU(5)®U(1) X 

/ (i.e., U(l)^(g>U(l) x -*U(l)g in the largo VEV limit.) 

(c) E 6 — > SU(3) C ® SU(2) L ® SU(2) R ® U(1) L ® U(1) R 

" v ' 

ER5M . , 

U(1) V=L+ R 

(d) E 6 — >SU(3) C ®SU(2) L ®U(1) Y ® SU(2) t ®U(1)' 

V v ' 

ER5M L SU (2) I 

Figure 1.6: E 6 Wilson-loop-breaking schemes jL3|. (a) shows a rank-5 model and (b) through 
(d) show rank-6 models. Scheme (a) gives the SM plus an extra U(1)y e - Schemes (b) through 
(d) can produce effective rank-5 models, ER5M, by taking a large VEV limit. 

The Supermatter Fields 

The most general superpotential that is invariant under SU(3) C (g> SU(2)l <2) U(1)y 
and renormalizable for the fields given in figure |1.5| is of the form (neglecting various 
isospin contractions and generational indices) |Tc 



W = Wo + Wi + W 2 + W 3 (1.2) 
Wi=A 6 $ dL $ u c$ e c + A 7 $ l $^$q + Ag^^^^c 

W / 3 = An$ R /$ L $ I/ //c . 

$ A = A) is the superfield, such that A = R', Q, ii£,. . . , and 



$o = I I , $r. = I I , $r, = I I , 
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for the first generation of the 27's, and similarly for the other generations. The 
Yukawa couplings, Aj's, also carry generational indices which have been suppressed; 
the couplings are inter-generational as well as intra-generational. The superpotential, 
W, summarizes the entire possible spectrum of low energy physics which can occur 
within the context of an Eq framework. 

Notice that W was only required to be invariant under the SM gauge group. 
Further constraints from E 6 model building may cause some of the Aj terms to 
disappear. Furthermore, not all of these terms can simultaneously exist without 
giving rise to AL^ and AB^ interactions; E 6 models say nothing about the 
assignments of baryon (B) and lepton (L) number until they are connected to SM 
representations. As a result various scenarios may occur, 

• Leptoquarks: B(q ! L )— | , L(q' L )—l =>- A 9 =A 10 =0 

• Diquarks : B(<£)=- § , L(<£)=0 =S> A 6 =A 7 =A 8 =0 

• Quarks : B(q' L )= ± , L(q' L )=0 =>- A 6 =A 7 =A 8 =A 9 =Aio=0 

where it has been assumed that L(z/f L )= —1 (these scenarios assume that there exist 
only three copies of the 27; more complicated ones can be constructed by adding 
extra copies). In this thesis the least exotic of these models, i.e., the "Quarks," will 
be investigated. Furthermore, to avoid any fine tuning problems with the neutrino 
masses, 

m v c « m e <t=^ An << A 3 , 

it will be assumed An=0. 

In this model the masses of the particles are generated by letting the role of the 
Higgs fields be played by 

L' = T eL , R/ = r L i ' r 



for each generation. It is possible to work in a basis where only the third generation 
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of Higgses acquire a VEV; the remainder become "unHiggses" [|T3j, [T^] . In this basis 
the Yukawa couplings, 

where k = 1, 2, 3 are generational indices, takes on a much simpler form, 

A 4 G {Af |Af 3 = Af 3 = A 33i = , A 333 = X 3 A jk = Af k = \f 3 ? s.t. i = 1, 2 & j, k = 1, 2, 3} . 

This basis also eliminates the potential problem of flavour changing neutral currents 
at the tree level. It is also assumed that the Aj's are real and that the couplings to 
the unHiggses are very small. The former assumption helps to further simplify the 
model and reduce any effects it might have in the CP violating sector JIB[ . 

1.2.2 Heavy Lepton Production 

E 6 models are very rich in their spectrum of possible low energy phenomenological 
predictions. If any new particles are found that fit within this framework then 
perhaps it will lead the way to a more unified theory of the fundamental forces of 
nature. However this is no small task, for a full theory would have to be able to 
actually predict the mass spectrum of the particles and the relationships between 
various couplings, and yet require very few parameters. Superstring inspired E 6 
models are far from being able to complete this task. However, proof that E$ is 
an effective GUT would be a good first step. But even this would not necessarily 
qualify superstrings to be the next step for it is not totally inconceivable that some 
other theory might give rise to E 6 as an effective GUT — caveat emptor. 

A natural question to ask would be, "Where to look for Eq phenomenology?" 
High energy hadron colliders, such as the Tevatron at Fermilab (1.8 TeV c.o.m., 
C ~ 10 2 p6 _1 /|/r, pp) or the LHC (14 TeV c.o.m., C ~ \Q b pb~~ l /yr , pp), offer pos- 
sibilities of observing phenomena beyond the SM by looking for the production of 
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q.q 




Figure 1.7: Gluon-gluon fusion to two heavy leptons, gg — ► L + L~ . The loop contains quarks, 
q, which couple to vector bosons, Z i= i i2 , scalar-Higgses, H° =12 3 , and a pseudo-scalar-Higgs, 
P° , and squarks, q , which couple to scalar-Higgses. 



heavy leptons through a mechanism known as gluon-gluon fusion, see figure |L7. 
This is an interesting process because there are enhancements in the cross-sections 
related to the heavy (s)fermions running around in the loop. The computation was 
done in the minimal supersymmetric standard model (MSSM) by Cieza Montalvo, 
et al., [16| in which they predict O(10 5 ) events /yr. Therefore for Eq it is expected 
that the production rate should in principle be higher since there are more particles 
running around in the loop. This process will be investigated in this thesis. 



1.3 Summary 

String-flip potential models and superstring-inspired-E 6 models have been discussed 
in a general setting. In this thesis various aspects of these models will be discussed. 
In chapter 2 the SU(3) C string-flip-potential model will be investigated and put into 
perspective with more generalized models. In chapter 3 modifications to the general 
string-flip-potential will be investigated. In chapter 4 heavy lepton production via 
gluon-gluon fusion will be investigated in an E 6 framework. Chapter 5 will contain 
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Chapter 2 

The SU(3) String-Flip Potential 
Model 

2.1 Introduction 

In this chapter a string-flip potential model for 3-quark systems will be constructed. 
Some simplifying assumptions about flux-tube minimization will be made in order to 
reduce the Monte Carlo computation time. The results for a linear potential model, 
SU^(3), and a harmonic oscillator potential model, SU/ l (3), in which the colour 
has been fixed to a given quark, will be presented. The results will be compared 
with an SU/ l (3) model proposed by Horowitz and Piekarewicz M in which different 
simplifying assumptions, about the minimal flux-tube topology, were made. Also a 
comparison of their results Q with some earlier work done by Watson || on SU(2) 
will be made. The chapter concludes with a discussion on possible future directions 
to pursue in order to obtain bound state nuclear matter. 
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2.2 The General String-Flip Potential Model 

As mentioned in the previous chapter a crude quark model of nuclear matter would 
be expected to have to at least the following properties: at low densities the quarks 
should condense out to form isolated baryons; at a higher density, the interaction be- 
tween quarks should lead to positive binding energy between nucleons and a swelling 
of nucleons; and at still higher densities, it is expected that the hadrons should dis- 
solve into a quark-gluon plasma. This last assumption is in contrast to the traditional 
nucleon models, which require the forces to be carefully adjusted so that they satu- 
rate at infinite density, effectively implying a hard core. Some simple models which 
appear to be likely candidates are string-flip potential models [Q, [5], and to some 
extent linked cluster expansion models [|l7j . 



The cluster models are based on one-gluon exchange potentials and use an N- 
body harmonic oscillator potential, i.e., 

Konf = \k £(?,- fjf, (2.1) 

to mimic quark confinement. These models are mainly used for describing short 
range nuclear effects, as they suffer from van der Waals forces due to the nature of 
the confining potential. Despite this shortcoming, they appear to be quite useful in 
explaining local effects such as nucleon swelling (fat nucleons) Jl8| , quark clustering 
preferences, and relative strengths of the various one-gluon exchange potentials fl7 



The string-flip potential models are, on the other hand, motivated by lattice 
QCD. These models postulate how flux-tubes should form amongst the quarks at 
zero temperature based on some input from lattice QCD. An adiabatic assumption 
is made, in which the quarks move slowly enough for their fields to reconfigure 
themselves, such that the overall potential energy is minimized: i.e., 

V = min{ v(r m ...r n )| |J {q m . . . q n } = {q x . . . q Nt }} , (2.2) 

{q m --qn} {m—n} 
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where the N q quarks are placed in a cube of side L and subjected to periodic bound- 
ary conditions to simulate continuous quark matter. The sum is over all gauge 
invariant sets {q m . . . q n } of quarks, such that at least one element from each set lies 
inside a common box, whose disjoint union, U , makes up the complete colour singlet 
set {gi . . . qN q } of N q quarks. It is easy to see that this potential allows for complete 
minimal quark clustering separability at low densities without suffering from van der 
Waals forces. At present, these models [f|, [| are quite crude in that they do not 
include short range one gluon exchange phenomena and spin effects and are flavour 
degenerate. Despite their shortcomings, in general, these models seem quite capable 
of correctly describing most of the bulk nuclear properties, with the exception of 
nuclear binding. 

It is known that the SU(2) string-flip potential models do show these properties, 
except for the positive binding energy which probably arises from short range forces. 
However, the only existing extension to an SU(3) model leads to the rather 
surprising result that the nucleon appears to shrink in nuclear matter. It is therefore 
of some interest to repeat the calculation of M, in an attempt to determine whether 
the approximations made there alter qualitatively the solution. 



2.3 SU(3) String-Flip Model 

The string-flip model involves solving a Hamiltonian system of fermions governed 
by the potential given in equation fl2.2|). To solve this system requires the use of 
variational Monte Carlo techniques |19 |. In order to compute any observable in a 



finite amount of time further assumptions about the form of the potential must be 
made. 
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In this model the potential is restricted 

a) to summing over sets of colour singlet clusters of three quarks, 

V = min{ v ( ? r,? g ,r b )| |J {q r q g q b } = {qi ... q Nq }} 

{irqgQb} {rgb} 

b) such that the colour of a given quark is fixed. 



(2.3) 



Assumption a does have some validity as it has been shown, via a linked quark 
cluster model, that it is energetically more favourable for Qq systems to dissociate 



into two nucleons as a result of hyperfine interactions [|17], ^TJ. However, this is 
not necessarily the case at lower densities, as the linked cluster models are unreliable 
here. Assumption b, that of fixed colour, greatly restricts the number of possible field 
configurations and therefore reduces the chance of finding an absolute minimum. At 
low densities this should not have any effect on the potential, as the system consists 
of isolated nucleons. Similarly at high densities no effect is expected, as the system 
consists of uncorrelated quarks. At intermediate densities some effects might be 
expected, particularly around any regions in which a phase transition might occur. 
For the SU^(3) model, the potential, V, has the components || 

r br + r rg if L brg > 120° 

r rg + r gb if L rgb > 120° 

r gb + r br if L gbr > 120° , (2.4) 



Wr,r r b j 



o < 



I V2 



ijse + vzA 



otherwise 



where r^ = r \ — Tj , £ or 



and 



£rgb — ^ \Z r ?g + r gb + r br > 



A — (r rg + r gb + r br )(— r rg + r gb + r br )(r rg — r gb + r br )(r rg + r gb — r br ) 



(2.5) 



(2.6) 
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is the area inclosed by the triangle Argb (see figure pTT|.a). For SU^) the compo- 
nents are 

v h (r r ,r g ,r b ) = ^k^ gh (2.7) 

(see figure [2.1| .b). This potential was obtained by replacing the linear segments of 
ve by springs when the quarks, which were assumed to be of equal mass, formed a 
triangle with interior angles less than 120° (cf. figure |2.5| .a). This analogue model 
is expected to have similar features to SU^(3) for s-wave (qqq) states. 



a) 



b) 




Figure 2.1: Flux-tube arrangements for the 3q cluster potentials a) vi and b) v^. 

The sum in equation ( |2.3| ) can be re-ordered by restricting it to run over all 
sets, {rgb}^, of quark triplets contained in a central box, such that the potential 
v(r r ,r g ,rb) is minimized with respect to all possible periodic permutations of the 
vectors {r r ,r g ,Fb}, with the constraint that at least one of the vectors lies inside the 
central box: i.e., 

V = mm{J2 w (rr,r g ,?b)} , (2.8) 

{rgb>L 
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where 



V{T T ,T g ,T h 



min{v(r r + Lk r ,f g + Lk g ,i h + Lk b )\ {k q ) a = -1,0, 1 
& at least one k q = 0} . 



(2.9) 



This means that for quark triplets a search of one box deep from the central box is 
required, giving a total of 27 2 possible permutations, in order to minimize a given 
v(r r , r g , r b ) (see figure P72"D . These permutations can be reduced to 3 by requiring that 
at least two sides of the triangle Argb , formed by a given permutation of quarks, 
be a minimum: i.e., 

v(r T , r g , r b ) = min{v(5r rg , 5r gh , r br ), v(r rg , 5r gh , Sr hl ), v(5r rg , r gb , 5r br )} , (2.10) 

where ^fy is the minimum distance vector between the points Fj and Tj, in a box of 
side L with periodic boundary conditions, which is given by 



(fi-T^a-L if (Ti-Tj) a > L/2 

where a — x,y, z . This is exact for SUft(3); for SU^(3), classical Monte Carlo shows 
that about 19% of the events deviate from the exact answer by ~ 0.3%, on average. 

The number of different elements in the set, {^}, from which the minimum 
must be extracted in order to obtain V, defined by equation (|2.8| ), is (iVJ) 2 (where 
N n = N q /3 is the number of nucleons in the central box). For example, if N n = 7 
this would yield 25 401 600 elements! These elements can be reduced by fragmenting 
the set {qi . . .qN q } into smaller pieces, or subclusters, such that each element can 
find Ntr complementary coloured pairs of quarks that are "closest" to it. These 
subclusters can be further fragmented, by "softening" the requirement that at least 



(^ij)a 



(ii-Tj) a + L if (fi-fj) a < -L/2 
(Ti-Tj) a if -?j)o| < L/2 



(2.11) 
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Figure 2.2: A 2-D slice showing a typical flux-tube arrangement for three quarks, placed 
inside a central cube, subjected to periodic boundary conditions. 
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Nbr complementary pairs exist: i.e., by searching for disjoint subclusters. These 
subclusters are referred to as softened subclusters. The "closeness" of quark q r to 
the complementary pair (q g qb) is defined by the function 5 r ,(gb) = v(f t , r g , Ft,), given in 
equation ( |2.9|) . The fragmented sets are thus constructed by computing an N n x N% 
matrix (A) with elements S r ^ g b), and then converting it into block diagonal form (A d ) 
increasing in size from top to bottom, by swapping rows and columns such that each 
block diagonal element contains the elements of a fragmented set and all the off 
block diagonal elements are set to zero. The elements of {X^} are now constructed 
by extracting permutations of elements S r ,(gb)> from unique columns and rows, of the 
block diagonal elements of A d . Further computational speed is gained by discarding 
sums that start to exceed the current minimum. In general, the fragmented sets 
constructed from {q% . . . qN }, are not all disjoint from one another, and therefore 
the block diagonal elements of A d may overlap. The degree of overlap increases with 
increasing density, causing the Monte Carlo to slow down. 

This fragmentation procedure, or nearest neighbour search of depth Nb r (cf. |jj), 
reduces computation time quite significantly. The cost is that rare configurations 
with flux-tubes that stretch across the box or across unsoftened subclusters, that give 
a global minimum, might be missed. Preliminary Monte Carlo results show that the 
inclusion of softened subclusters yields no noticeable change. However, for a full 
(N n \) 2 brute force search, performing a Monte Carlo becomes virtually impossible. 
A few brute force computations of the potential were made, for particles randomly 
thrown into a box, which seem to suggest that the fragmentation procedure is good 
to about 1% , with Nb r ~ 4. SU^(2) models also give similar results ||. 

The validity of the fragmentation procedure can also be argued on physical 
grounds, for it is reasonable to assume that long flux-tube configurations would 
tend to dissociate into qq pairs. Therefore, the fragmentation procedure can be 
considered as a valid low temperature approximation. 
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The choice of variational wave function should attempt to reflect the overall bulk 
properties of the system. Here the wave function was chosen to be of the form, 

- e cur 

*(a,/3,p) = e ^b} Y[ |<Mp)| , (2.12) 

ce {rgb} 

where a, /3, and p (= N n /L 3 s.t. N n = N g /3) are variational parameters, Z){ rg b} is 
over the set of quarks {rgb} which gives the minimal potential V, and |$s c (p)| is 
the Slater determinant which is a function of density, p.0 The $s c (p) contain the 
elements faj = <fii(?j), which are composed of the plane wave states 

27T 

sin(— Ui ■ fj + 5i) , (2.13) 



where Si = or 7r/2 , and (hj) a = 0, ±1, ±2, . . . are the components of the Fermi 
energy level packing vector, hj, for particles in a cube of side L (with ordinates rang- 
ing from —L/2 to L/2) subjected to periodic boundary conditions. The exponent 
part to the left Slater part of equation |2.12| yields a product of three-body harmonic 
oscillator wave functions when a = 2. Therefore, the parameter (3 is related to the 
inverse r.m.s. radius of the nucleons in the system (cf. ||). Our particular choice 



of variational wave function, i.e., equation |2.12| , mimics the overall gross features of 
quark matter, yielding highly correlated behaviour at low densities and uncorrelated 
behaviour at high densities. 

The total energy for this many body system is, 

E(a,(3,p) = T^ + V } (2.14) 

where 

T_ s = ^(V 2 lnvI/) (2.15) 
Am q 



§ B.l gives an algorithm for generating an arbitrary dimensional Slater determinant. 
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is the kinetic energy, obtained by eliminating the surface terms from the integral 
J \]>*v 2 \I/ . Thus, the many body Hamiltonian system can be solved by varying 
the parameters a, (3, and p, and evaluating the expectation values by Monte Carlo 
integration at each step, until a minimum E is found. 



2.4 Monte Carlo Calculation 



The Monte Carlo procedure uses the Metropolis algorithm [IS, 22, 23, 24} to generate 
a distribution in |\l/| 2 . The Monte Carlo procedure involves computing the average 
of an observable O such that 

1 r i N+N 

6 = j^r- j 0(x*)|*(x)| 2 dx « - £ 0(x n ) , (2.16) 

where the summation is taken over N sequential samples of the distribution \^(a, 
j3, p;x)| 2 (s.t. a, j3, and p are fixed), after N iterations have been made. The dis- 
tribution in |\I/| 2 is generated by the following algorithm: from a given configuration 
of particles x, change all their positions randomly to a new position x + dx . Next 
compute the transition probability function 

t = min{(|^(x + dx)|/|^(x)|) 2 ,l} , (2.17) 

and compare it with a random number re [0, 1) . If r > r then accept the move by 
replacing x with x + dx, otherwise reject the move by keeping the old x. Finally, 
repeat the procedure until a desired error level, 50 , has been reached. The initial 
configuration of particles, x = x , is generated by throwing them randomly into a 
box with side of length L. All subsequent moves are constrained to the box such 
that, if a particle randomly moves outside of the box, its periodic image enters from 
the opposite side. This algorithm, which satisfies detailed balance, is called the 
Metropolis algorithm and converges to the distribution |\1/| 2 after iVo moves have 
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been made. The value of Nq is determined by the point at which the statistical 
fluctuations in J2 n @n have become substantially reduced. A general rule of thumb 
is that convergence is more rapidly achieved if the step size, 5 s x (= \dx\/y/3N q ), 
is chosen such that, on average, r ~ 1/2. A natural length scale to use, when 
considering an appropriate step size, is (cf. f|) 

where the constant / ~ 1/4. This is determined by taking several small samples 
from the probability distribution |\l/| 2 and by restarting the Monte Carlo for different 
/ values, until the desired value of f is reached. To ensure convergence in a finite 
amount of cpu time, particularly at low densities, rgb clusters of quarks (of radius 
order 5 s x/2) are thrown into the box randomly. 



The Monte Carlo evaluation of the total energy, E = T_ s + V in its current form, 



can produce a significant amount of error [IS]. This can be reduced by introducing 
a mean square "pseudoforce" , 

F2 = ii™ 2 - < 2 - 19 > 

and re-expressing the kinetic energy as 

T = 2T„ S - F 2 . (2.20) 
In this form, the variance of the total energy, 



E = 2T„ S -F 2 + V } (2.21) 



goes to zero as the wave function approaches an eigenstate of the Hamiltonian. 

The variational wave function \& is made up of a product of a correlation piece, 
X, and a Slater piece, $. Therefore, the kinetic energy expression can be split up 
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into three separate terms involving pure and mixed, correlation and Fermi energies: 
i.e., 

f = f c + T F + f CF (2.22) 
The explicit forms for these terms are: the correlation energy 

nP> a 

To = ^-V< E Cb 2 Hl - (/^gbD + 4]) , (2.23) 

Z?7 V V ™ {rgb} 



the Fermi energy 



27T 2 V 



f F = — -— E 5 '' (2-24) 
m q N n L l p[ 



and the mixed correlation-Fermi energy 

T ° F = E 4gb [0ir( ? rg - ?br)0- r + 0ig(?gb - ? rg )0- g 

i=l {lgh} 

+0ib(?br - ? g b)0'ib] • ni) , (2.25) 

where </>y = {4> T )ij 1 , and 0'^ = 0^ + ^rij) . A detailed derivation of the above 
expressions can be found in appendix |A]. 

The value of a is fixed for free nucleons at p = 0. For the SU^(3) model a = 2, 
as the wave function \1/ must become that of a free 3-body harmonic oscillator. 
Therefore, the total energy for this system is simply 

rf&W = ?V + H • (2-26) 



Minimizing this gives, 



E% 1 > = 3h 



\ 



_k_ 

m q ' 



(2.27) 



where E^ = E^ e (p { h) ), and 



0g>= raj • (2-28) 



1/4 
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E [ h) and (3 { h) can be used to check the Monte Carlo. However, for the SU£(3) model 
such a check is not possible, as it is impossible to find V analytically at p = . 
However, by fitting the results to the expression 

*%LW=!h— P 2 + 9v^, (2-29) 

we find g T ~ 1.07 and g v « 3.09. Also the virial relation (7^) = (Ve)/2 can be 
verified, which should hold at all densities. A similar check can also be done for 
SUft(3), with the virial relation (Th) = (Vh) ■ For SU^(3) the parameters a and 
can be obtained by computing E^l e for different values of (a,/3) until a minimum is 
found. 

A further reduction of the variational parameters is obtained by introducing the 
scaling transformation H, 

(/3,p 1/3 ) — > C(cos0,sin^), (2.30) 

where ( > 0, and 9 is restricted to the interval (0, 7r/2). This allows the total 
energy to be expressed as a polynomial in (, which can subsequently be minimized 
to eliminate (: i.e., E(p,/3) becomes 

E((,9) = f(e)( 2 + V(9)/C, (2.31) 
such that f (C, 9)=f (9)C , V"(C, 9) = V{9)/( K and 

f 1 ifSU/3) 

j (2.32) 

[ 2 if SU h (3) 

which can be minimized with respect to £ to give 

gW = (« + 2)|' ^" W | . (2-33) 
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with 



kV(9)\ k 



< + 2 



(2.34) 



2T (9) J 

Notice that the elimination of the parameter ( is equivalent to imposing the virial 
theorem, which implies 

(T) = ^(V). (2.35) 

Therefore, the Monte Carlo only has to be run for different 9 values extracted 
from the "open" interval (0,7r/2). The end points are obtained by taking a limit. 
The 9 = limit is equivalent to taking p = , which has already been discussed. The 
9 = 7r/2 limit is equivalent to taking (3 = 0, which corresponds to an uncorrelated 
Fermi gas, with energy 



E Fg ( P ) 



5^4 \ 1/3 



3 5 7T 



3tr 

5m, 



p 2/3 + V Fg (p) 



(2.36) 



where 



V Fs ( P ) 



C K 



a 

^73 



k 



for SU £ (3) 



for SU fe (3) 



(2.37) 



and c K is obtained by a fit to the Monte Carlo in the 9 = n/2 limit. Thus the (3 = 
limit is described by the curve E Fg (p) . This curve is compared with the Monte Carlo 
results for E(p(9)) from which a minimum energy curve E(p) is obtained. 

Figures ^]3| and ^]4| show the variational Monte Carlo results for (3(p) and the 
binding energy, Eb{p) = E(p) — E , respectively. The dashed lines on these graphs 
show the remnants of the minimal p(9) trajectories, for (3 and Eb , after a phase 
transition, at p = p c , from a correlated system of quarks to an uncorrelated Fermi 
gas. The slight roughness of these lines is because the data was not fitted. In plotting 
these graphs it was assumed that: N n = 7 , m q = 330MeV , a = 910MeV/ fm , and 
k m 3244MeV/ fm 2 . The value of k was determined by setting E^ in equation 
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Figure 2.3: Graph of (3{p) for SU^(3) and SU fe (3) . 
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Figure 2.4: Graph of E B {p) for SU £ (3) and SUfc(3) . 
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Table 2.1: Parameters determined by Monte Carlo, with N n = 7 . 



Parameters 


PITT /O \ 

bU*(3) 


PITT / 0\ 


a 




1.75 


2.00 


c k(Pf s , Pf s ) 




0.908 


0.347 


Pf s 


(fm 3 ) 


0.2524 


0.465 


f3 Fg 


(MeV) 


0.01264 


0.01549 


E (p ,p ) 


(MeV) 


1895 


1856 


Po 


(fm- 3 ) 


1.119 x HT 11 


1.211 x 10~ n 


Po 


(fm- 1 ) 


2.237 


2.297 



( |2.27| ) to equal E$ , in the limit = 1 and = 0.0001 . The other parameters 



that were determined by the Monte Carlo are given in table ( 2.1 



2.5 Discussion 



The parameter (3 is related to the confinement scale for triplets of quarks ||, |5|. 
Figure [T3] shows that the quarks become less confined as p increases, and completely 
deconfined beyond the phase transition point, p c . Thus, as p increase from to 
p c the nucleon swells producing an EMC-like effect. Figure 11. b, of reference 0, 
shows a plot of A (~ f3 2 ) vs. p, obtained by Horowitz and Piekarewicz, which in 
general indicates that the quarks become more confined as p approaches p c , and 
completely deconfined beyond. Therefore, their model does not explain the EMC 
effect [E25|, 



, |2]J and is inconsistent with what we have found here. 
The Horowitz and Piekarewicz (HP) model approximates the higher order flux- 
tube topologies of equation (|2.2|) with long harmonic oscillator chains that close upon 
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themselves: i.e., 



V = min{^ f rg } + min{^ t> gb } + min{^ v hr ] , (2.38) 

{rg} {gb} {br} 



where Vij = \krf-, with the wave function \I/ = e~ Ay <l> (cf. equation ( 2.12 ) with 
a = 2). They have shown for p < p c that 3-quark clusters (chains) make up more 
than 90% of nuclear matter, with a large remainder of these being 6-quark clusters. 
A closer look at their A vs. p plot shows a small dip in A around p = 0.2 . This 
indicates a slight swelling of the nucleon for very small p. In fact, in this density 
regime 3-quark clusters completely dominate (> 99%). Therefore, the evidence from 
these graphs seems to suggest that too much weight is being given to higher order 
flux-tube topologies at intermediate densities. 

Lattice QCD shows that quarks like to cluster together via a linear potential. 
However, most phenomenological models that describe isolated hadronic matter us- 
ing a harmonic oscillator potential work just as well. As can been seen from figures 
273] and [Z~4] the harmonic oscillator model gives the same overall shape as the lin- 
ear one. This model was motivated by replacing each linear segment of string in a 
3-quark state by a spring, with spring constant k. For quarks of equal mass this re- 
duces to a triangle of springs (see figure p75[a). Similarly a 6-quark state would give 
an object that simplifies to three triangles with one of the tips from each meeting 
at a common vertex (see figure |27J.b). The corresponding 6-quark state for the HP 
model forms a closed ring which, in general, requires less energy to form (see figure 
p.5| .c). Thus QCD motivated models would also seem to support the aforementioned 
claim, that HP are giving too much weight to higher order flux-tube topologies. 

The SU/!(2) model, by Watson 0, agrees with the HP model ||]. These graphs 
appear to be similar to those shown in figures |27^ and [2.4| . Of course the fact that 
these models agree is only a check of consistency, as they both have the same poten- 
tial, which only includes interactions between qq pairs. Also the SU/i(2) models |4], [5[] 
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Linear Harmonic Osc. 




Figure 2.5: Different SUh(3) flux-tube construction schemes (RHS) motivated by their cor- 
responding linear cousins (LHS). Figures a) and c) represented the HP construction. Figures 
a) and b) show a construction scheme with a more consistent weighting for s-states. 
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when compared with SU^(2) || gives similar contrasting figures to the ones presented 
here. 

Figure 2.4 , along with similar figures given in references ||, show a saturation 



of nuclear forces as p — > p c , followed by a phase transition to quark matter at p c . 
All of these models, however, fail to give any nuclear binding below p c , which would 
seem to suggest that the flux-tube models are incapable of obtaining nuclear binding. 
Even the SU^(3) HP model with its long chains, which tends to underestimate the 
potential, indicates that this would appear to be the case ||. HP have a 2q model 
H that would seem to suggest that even if colour were not fixed to a given quark, 
no nuclear binding would occur: albeit this model is for p-wave (qq) states. Thus, it 
would appear that string-flip models, even those that include higher order flux-tube 
topologies, or allow the colour to move from quark to quark, are insufficient to obtain 
nuclear binding. Therefore, another mechanism for obtaining nuclear binding must 
be included in these models. 

One possibility is to include one-gluon exchange interactions. As suggested by 
Nzar and Hoodbhoy ||17|| , the most significant of these are the hyperfine interactions. 



In a relativistic setting one could also consider the possible effect of chiral symmetry 



breaking in which the constituent quark mass changes with momentum scale p8| , [29|j.a 
Finally, other effects such as quark mass differences and isospin are expected to be 
negligible. 

2.6 Conclusion 

Various string-flip potential models have been discussed in general, and have been 
shown to adequately describe the bulk properties of nuclear/quark matter with the 



2 This interesting possibility was pointed out to us by the referee of reference 
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exception of nuclear binding. At low densities they yield free nucleon matter and 
at high densities a phase transition to free quark matter. They show an overall 
saturation of nuclear forces as nucleon densities are increased. At intermediate den- 
sities these models, with the exception of the HP SU/ l (3) linked chain model [[|, give 
an overall EMC-like swelling of the nucleon. It is believed that these models with 
the addition of one-gluon exchange effects should be capable of predicting nuclear 
binding. This is the topic of the next chapter. 



Chapter 3 

Flux-Bubble Models and Mesonic 
Molecules 

3.1 Introduction 

In the previous chapter string-flip potential models were investigated. It was hypoth- 
esized that these models by themselves were incapable of producing nuclear binding 
and would therefore have to be extended. The suggested extension was to include 
one-gluon exchange interactions. In this chapter a new class of models, called flux- 
bubble models, is proposed which allows for the extension of the flux-tube model to 
include these interactions. 

3.2 Flux-Bubble Models 

The primary objective is to construct a model which combines both nonperturbative 
(flux-tubes) and perturbative (one-gluon exchange) aspects of QCD in a consistent 
fashion. In order to simplify this task only the colour Coulomb extensions to an 
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SUe(2) potential model will be considered. 

The extension of the linear potential model for qq pairs is 
cr(r.y - r ) if T i:j > T v 



aSij [^~7 Q ) lfTij<T ° /^Coulomb 

where Ajj = —3/4, 1/4 for unlike and like colours respectively, and a s ~ 0.1 . This is 

simply a variant of the phenomenological potential, 

4 rv 
3 r 

mentioned in chapter 1; the major difference being that the nonperturbative and 
perturbative parts are completely isolated in the former as opposed to the latter. 
When the quarks are separated at a distance greater than r the potential is purely 
linear and when they are inside this radius it is purely Coulomb. In effect, for 
distances less than r , a "bubble" is formed in which the quarks are free to move 
around, in an asymptotically free fashion. In both distance regimes the net colour 
of the system is neutral. 

The extension of the linear potential although simple for a pair of quarks becomes 
more complex when considering extensions for many pairs of quarks. In particular, 
how can a potential model be constructed in which some quarks are close enough 
to be inside perturbative bubbles while at the same time a subset of them are still 
connected to nonperturbative flux-tubes that extend outside of these bubbles. The 
solution to this problem is easily remedied by inserting virtual qq pairs across any 
of the intersection boundaries formed by the flux-tubes with the bubbles. Now the 
segments of flux-tubes that lie outside the bubbles remain intact while the segments 



inside simply dissolve; giving the desired result. Figure |3.2| illustrates the dynamics 
of this model. Notice that this model allows the construction of colourless objects 
because of the insertion of the virtual qq pairs. These virtual quarks are used as 
a tool to calculate the overall length of the flux-tube correctly. They are not used 



CHAPTER 3. FLUX-BUBBLE MODELS AND MESONIC MOLECULES 



41 



Figure 3,1: Consider configuration (a) of quarks, with r > ro, 

a) m q 

about to move to (b), s.t. two of them are within r < ro. Then ^ 
the procedure is to draw a bubble of r away from the two, (b), 
and to cut the flux-tubes at the boundary and insert virtual qq 
pairs, (c). Once the potential is computed the configuration is 



restored to (b) before the next move is made. 

b) c) — Nonperturbative 

Piece 






Cut here 



Virtual qq pair insertion 



in computing the Coulomb term however, as the field energy is already taken into 
account by the "real" quarks inside the bubbles. In general, once the bubbles have 
been determined, the flux-tubes must be reconfigured in order to minimize the linear 
part of the potential. 

Although the model is currently for SUi(2) it should be easy to extend it to a 
full SUe(3) model with all the one-gluon exchange phenomena. 



3.3 In Search of a Wave Function 



Figure |3]2| shows some preliminary results using the SU^{2) flux-bubble model de- 
scribed in the previous section. For ro = fm the Monte Carlo recovers the same 
results for the SU^{2) string-flip potential model ||, as expected, and for r = 0.1 fm 
the result differs only slightly. 

These results are questionable, as the wave function that was used is not ideal. 
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Figure 3.2: E B (p), (a), and /3(p), (b), for r = Ofm 



and rp = 0.1 fm( — ), with 



m q = 330MeV , a = 910MeV/ fm , and a s = 0.1. These graphs were created by using a 
coarse 10 x 10 mesh of points in p and p. The minimal curves were extracted by using linear 
interpolation between the minimum data points on the Eb (p, /3)-mesh surface. 

It consisted of a slight modification to an old SUg(2) wave function ||, 

- E (0r„ ^ 



aa'fip 



min{qq} 



(Fermi) 



(Linear) 
t a 



cf. equation |2.12|) , in which a new correlation piece, x a ' 



(Coulomb) 



(3.2) 



was added to account 



for the local attractive Coulomb interactions as they occurred: i.e. 



j (Linear) (Coulomb) 

Vaa'Pp ~ X Q /3 X X a 'f3 



(Fermi) 
P 



(3.3) 



where 



(Coulomb) _ min{qq} 
Aa'/3 ~ e 



(3.4) 



with a = 1.75 and a' = 1 (the repulsive interactions were assumed to be taken care 
of by the presence of the Slater wave function, $( ) Fermi ) ) . However, because the flux- 
tubes and bubbles can now be created or destroyed the wave function is, in general, 
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no longer continuous to order V 2 (cf. equation |3.4| ), and so a variational lower bound 
is no longer guaranteed. Therefore, a new wave function is needed. Unfortunately, 
it is a rather difficult task to come up with a wave function that takes into account 
the locality of these flux-bubble interactions which is smooth and continuous, and 
involves very few parameters. 

The aforementioned wave function has two "independent" parameters, p and 
(3; the a's were assumed to be fixed. This is in contrast to the previous case for 
the SUe(3) model, in chapter |2|, in which p and (3 varied parametrically with a 
single parameter, 9 . The reason this does not apply here is because the flux-bubble 
potential breaks the scaling transformation 

(Ap 1 / 3 ) — C(*)(cos0,Bin0). O) 
This extra degree of freedom greatly increases the computation time. 



The results in figure |3.2| were generated on a coarse 10 x 10 mesh of points, in p 
and P, and required 18hrs of CPU time on an 8 node farm. Clearly this procedure 
would be ridiculously slow if more parameters were to be added: 

r CPU ~ 0(R P ) (3.5) 

where R is the mesh resolution in each ordinate, and p is the number of parameters. 

Therefore, a way of checking different wave functions and minimization schemes 
which does not consume large amounts of CPU time is desirable. In particular, a 
mini-laboratory is needed in which various aspects of the string-flip and flux-bubble 
potential models, from wave functions to minimization schemes, can be investigated. 



3.4 Mesonic Molecules 



Some work was done with mesonic-molecules by Treurniet and Watson [31], 
using SUe(2), which has shown that these molecules make useful mini-laboratories for 
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studying string-flip potential models. They used a mesonic-molecule, Q2, consisting 
of two heavy quarks and two relatively light antiquarks: see figure |3l| The quarks 




q 



Figure 3.3: A mesonic-molecule, Q2 , with two heavy quarks and two light anti- 
quarks. 

are assumed to be heavy so that the light antiquarks can move around freely without 
disturbing their positions. By varying the distance, R, between the heavy quarks 
a mesonic-molecular potential, U(R) , can be computed (cf. equation |3.11| ). The 
Q2 system provides a good way of checking potential models for the possibility of 
nuclear (mesonic) binding. Moreover, because of its simplicity, it allows for checking 
various wave functions and minimization schemes without being concerned about 
CPU overhead. 

3.4.1 The Distributed Minimization Algorithm 

There are many ways of determining the minimum of a function, f(x), of several 
variables, x — (xi,. . . ,x n ) [0. However, for the case where the function is approxi- 
mated by Monte Carlo most of these methods, in general, will not work. The reason 
for this is because Monte Carlo calculations produce results that have statistical 
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uncertainties; most of the methods are for functions that give "exact" answers. For 
example, if gradient methods were used, then the error would propagate into the 
gradient calculations which would effectively add noise to the search. 

Of the viable methods the most promising one is the parabolic minimization 
algorithm f33|. The algorithm is as follows (see figure p.4|) : 



• Pick an initial starting point, Xj . 

• Pepper its neighborhood with p points, x p . 

• Evaluate the d — D function, f(x), at these points. 

• Fit a d — D parabola to these evaluated points, (f(x p ,i), x Pi i) . 

• Find the critical point, x c , of the parabola. 

• Use this point to repeat the procedure. 

This algorithm will, in general, converge to a critical point on the surface that is 
being searched. However, this does not guarantee that the point will be a minimum, 
it could very well be a maximum or a saddle point. Therefore, it is important to 



check the curvature |33| and then take appropriate measures to drive the search away 
from this point if is indeed not a minimum. Additionally, it is probably a good idea 
to bound the search to a box so that it does not drift out to infinity. It is important 
to note that this algorithm does not guarantee convergence to a global minimum but 
then neither does any other algorithm. 

This algorithm must now be adapted to take advantage of a distributed comput- 
ing environment. 

The first step is simply to submit p points to m computers; where if p < m use p 
computers and if p > m distribute p/m jobs per computer (a slight improvement on 
this last step would be to distribute m points and then dole out the remaining m + i 
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Local Parabolic 




Minimum. 



Figure 3.4: The parabolic minimization algorithm for an arbitrary (i-dimcnsional, d — D , 
energy surface, where p = \d(d — 1) + 2d + 1 is the minimum number of parameters required 
for fitting locally to a d — D parabola. 
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(i=l y . . ,p — m) points to the machines that finish first). However, this method is 
very inefficient because each iteration runs in time 



t ~ Q{ 



p_ 

rn 



(3.6) 



A much more efficient method would be to develop a procedure that analyzes the 
data in a continuous fashion as it streams in, point by point. 

Such an algorithm was developed over the course of time that it took to produce 
the results in the later sections of this chapter. The details of the algorithm can be 



found in § B.2. A basic outline of the procedure is as follows: 



1. Create an initial data base of k — m sample points, (f(x K - m ), x K - m ) , with m 
points pending. 

2. Fit the neighborhood of each sample point to a parabola. 

3. Submit the critical points, x c , that correspond to parabolic minima. 

4. Keep the m computers occupied by submitting extra sampling points if neces- 
sary. 

5. Wait for a point, (f(x), x) , to arrive. 

6. Update the data base. 

7. Fit the point about its neighborhood to a parabola. 

8. Submit the newly predicted point, x c , if it corresponds to a parabolic mini- 
mum. 

9. Go back to step 4. 
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This algorithm is effectively performing several parabolic minimizations all on dif- 
ferent regions of the d — D surface f(x) by interlacing its searches. For the work 
done here the convergence turned out to be quite rapid (due to cross-talking) with 
a fairly small start up cost, 



K 

m 



(3.7) 



where typically 

K ~ 0(2p)+m : (3.8) 

i.e., there were m points pending at the time of creation of the data base in step 1, 
above. Also the time between iterations (i.e., "events") was observed to be 

7*~0(-), (3.9) 
m 

where m < 8 . Ultimately this algorithm will saturate when approaches the time 
required to process each calculation, r p . r p is not a constant but a function of the 
size of the ever building data base and therefore will grow with time. Regardless 



though t p « Ti and therefore for all practical purposes equation [J]9] holds. At 
some point m will become large enough that this algorithm would begin to have the 
appearance of simple grid search. In fact, it was designed with this in mind. 

For the Monte Carlo results contained herein the number of iterations required 
to obtain an accuracy well below the 1% level was about 30, 60, and 150 for the 1, 
2, and 3 dimensional searches, respectively. Upon further investigation it was found 
that the convergence criterion that was used was too weak. The minimum was found 
after about 65% of the iterations needed to meet the convergence criterion. It is im- 
portant to point out that this algorithm is still in its infancy and requires further 
development. Indeed, later work with the 1 — D searches have reduced the itera- 
tions by about 50% and is expected to do the same for higher dimensional searches. 
The improvements were mainly due to establishing good convergence criterion and 
weighting techniques (see § |B.2| ) . 
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3.4.2 A General Survey of Extensions to SUg(2) 

In this section the effects of extending the old SU(_{2) || model to include flux- 
bubbles, with and without fixed colour, will be investigated in the context of the 
mesonic-molecular system, Q 2 ■ To simplify the situation these extensions will be 
considered in the frame work of SU(2) colour, SU C (2). 

A good place to start is by studying a Q2 effective potential, U (R) , that was 
generated by interactions with the light antiquarks through a linear potential, 

V = a r Q^ ( 3 - 10 ) 

min{Qq} 

where the sum, V] , is over the set of quark- ant iquark pairs, {Qq} , which min- 

min{Qq} 

imizes the potential, tq^ represents the distance between a given light antiquark, 
q , and a given - fixed - heavy quark, Q , and a w 910MeV/ fm . The Schrodinger 
equation that describes the effective potential between the two heavy quarks, in the 
adiabatic approximation, is as follows |34| 



2 ^v;- + y|* = w, (3.11) 

with m q ~ 330MeV . For the old SUe(2) model the variational wave function was 
assumed to be of the form 0, 

^i ap = e mm{Q 5 } . (3.12) 



Therefore, following a similar procedure to that of section |2.4| the effective potential 
is found by evaluating 

U = f c + V (3.13) 

where 

T c = a(3<*±-{ £ [ a (l-((3r Qq r) + l]r a Q f), (3.14) 

1 min{Q</} 
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(cf. equation |2.23| ), at different values of R. 

Figure [3.5| shows a plot of the Monte Carlo results for U(R) where a and (3 
have both been allowed to vary. This potential has been parameterized by a 
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Figure 3.5: The Qi potential obtained using a linear interaction with a pseudo- hydrogen 



wave function, equation 3.12 



pseudo-Morse potential of the form 

U(r) = { e - 2 ^-nO] Q _ 2 V e~ mr ~ ro)]a } a (r - r' ) + . (3.15) 
The term in the braces, "{}" , on the left along with a is the Morse potential p4| , [35 
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with two extra parameters included; r] and a. Outside the brackets to the right is a 
linear term with an offset parameter r' Q (this mainly takes effect in cases where the Q2 
potential blows up at the origin). Finally is the total energy at infinite separation. 
For these figures and the figures which follow, it will be assumed that rj = 1 and 
r' Q = unless specified otherwise. A summary of the values of these parameters can 
be found in table |3.1| . This parameterization was used mainly because it gave a fairly 
universal description of the (fully minimized) potentials in this chapter, however its 
physical significance should be taken lightly. 



Table 3.1: Summary of parameters for the fits in figure 3.5 



Potential 


Parameters 


Pseudo- Morse 


(3 (fm- 1 ) 


r (fm) 


a 


1.162 ±0.083 


-0.41 ±0.18 


3.0 ± 1.2 


a (MeV/fm) 


Uoo (MeV) 


r' (fm) 


11.9 ±7.1 


1709.4 ±0.4 


0.033 ± 0.060 


Yukawa 


V (MeV) 


a (fm) 


Voo (MeV) 


39.8 ±3.2 


0.3092 ±0.0075 


1709.40 ±0.02 


Parabolic 


C (MeV/fm 2 ) 


r (fm) 


Vo (MeV) 


32.6 ±2.8 


0.4060 ±0.015 


1706.2 ±0.1 



The values of Monte Carlo results for U(R) at the end points of the curve, from 
R = Ofm and out to R = 5fm, were checked against the analytic solution (cf. || 
and equation |2.29|) 

E free = 2(^(3> + g L (a)^) , (3.16) 

where 



9t((*) 



2-2 
a 2 2 a r(2 + l/a) 

r(3/a) 



(3.17) 



CHAPTER 3. FLUX-BUBBLE MODELS AND MESONIC MOLECULES 



52 



9l(<x) 



r(4/a) 



(3.18) 



2 a r(3/a) 

and /i (« m g ) is the reduced mass. This solution was easy to obtain because at the 
two extremes, R = and R = oo , V becomes separable; the Q2 system becomes 
equivalent to two isolated mesons. Equation |3.16| can be minimized in a and (3 easily 
by using Mathematica [35| . The results can be found in table [TP] of § |3.4.1| . 

It is interesting to investigate whether or not the string-flip potential, equa- 



tion |3.10| , actually mimics pion exchange. Therefore, the asymptotic part of U (R) , 
in figure |3.5| , has been fitted to a Yukawa potential: 

e -r/a 



V(r) = -V 



r/a 



+ V ol 



(3.19) 



and so the mass of the exchange particle is 



he 

a 



(3.20) 



A summary of the values of the parameters Vq , a , and Vqo can be found in table |3~I. 
Therefore, 

m cx w (638 ± !5)MeV , (3.21) 

which is about 4.6 times too big. 

Finally a test of whether or not the system is capable of binding can be done by 
expanding about the minimum of U(R) by using the parabolic approximation |3"5 



y (r) = C(r - r ) 2 + y 



which can be transformed into the harmonic oscillator potential 



V h (r) 



-D+ \li Q uj 2 (r - r ) 2 , 



(3.22) 



(3.23) 



with 



D = U oo -y , 



(3.24) 
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uj = c y2C/ fiQ , and fiQ = m,Q/2 , where Uoo shall be taken as the analytic result for 
U(R) out at oo. Therefore, the binding energy is 



E% = -D + hu;(v+l), 



which implies 



HQ 



> 



C(hcf 



2(U oo -y y 



(3.25) 



(3.26) 



in order to obtain binding. A summary of the parameters C , r , and yo , and the 
analytic result can be found in tables [D] and [31], respectively. So to obtain 
binding 

A*q £ (55 ± 6)GeV ; (3.27) 
which is not surprising as the well is quite shallow, 



D ss (3A±0A)MeV . 



(3.28) 



The next step is to add an SU C (2) term to the linear potential, equation |3.10 , 

(3.29) 



v= E 

min{Qq} 



o (r Qq - r ) 6(r Qq - r ) - ^ a s f -i- - ^-j #(r - r Q5 ) 



where a s ~ 0.1 , and tq « 0.1 /m . Figure |3]6| shows the results of the Monte Carlo 
for this potential. 

The Monte Carlo results for U(R) were checked at R = Ofm and R = 5/m 
against the analytic solution 

'#z(a,/3) 



free 



2[i 



g (a,f3)r 



3 ( , a* (l-9o(a,P)y 
•- a s U c (a, /3) /3 - 



(3.30) 
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Q 2 Potential Obtained by Using exp(-(/Sr)°) with V qq (r) = (lin) + (Coul). 
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Figure 3.6: The Qi potential obtained using a linear-plus-coulomb interaction with 
pseudo-hydrogen wave function. 
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where g T (ct) is given by equation p. 17 

r(4/a,2(/3r ) Q ) 



2 a T(3/a) 
T(2/a)-r(2/a,2(Pr ) a ) 

2 a T(3/a) 
T(3/a,2((3r ) a ) 



(l-P(4/a,2(/3r ) Q ))(? L (a), (3.31) 
P(2/a,2((3r ) a )g c (a), (3.32) 



r(3/a) 

where g L (a) is given by equation |3.18 



l-P(3/«,2(/?r ) ( 



(3.33) 



9c( a ) 



T(2/a) 



(3.34) 



2 a T(3/a) 

P(a, 2) = 1 — r(a, z)/T(a) , and T(a, z) is the incomplete gamma function [3B, 
At i? = 0/m and R — 5(~ 00) the system effectively decouples into two isolated 
mesons, for which equation [3.30| was derived. Notice that in the limits tq — > and 
r — > 00 equation |3.30| simplifies to the earlier linear solution, equation |3.16j , and to 
the Coulomb solution, 



Et — 2 



^/3 2 -|a s9o (a)/3 N 



(3.35) 



respectively. The analytic (i.e., the minimum of equation |3.30| , via Mathematica 
verses the Monte Carlo results are summarized in table |3~^ . 



Table 3.2: Monte Carlo (MC) vs. analytic results for a linear-plus-coulomb potential. 



Parameters 


Analytic 


MC @ R = Ofm 


MC @ R = 5fm 


E min (a,(3) (MeV) 


1527.07 


1527.58 ±0.04 


1527.61 ±0.05 


a 


1.74 


1.74 


1.72 


(3 (/m- 1 ) 


1.37 


1.37 


1.38 



The results for the pseudo-Morse, Yukawa, and parabolic fits are summarized in 
table O. 
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Table 3.3: Summary of parameters for the fits in figure 3.6 



Potential 


Parameters 


Pseudo-Morse 


P (fm- 1 ) 


r (fm) 


a 


2.01 ± 1.3 


-0.01 ±0.41 


2.0 ± 1.3 


a (MeV/fm) 


Uoo (MeV) 


r' (fm) 


11.2 ±6.0 


1527.6 ±0.5 


0.027 ±0.072 


Yukawa 


V (MeV) 


a (fm) 


Voo (MeV) 


13.61 ±0.55 


0.4299 ±0.0082 


1527.70 ±0.02 


Parabolic 


C (MeV/fm 2 ) 


r (fm) 


Vo (MeV) 


33.1 ± 1.6 


0.419 ±0.053 


1524.2 ±0.1 



The Yukawa potential gives an exchange particle mass of 

m cx ss (459 ± A)MeV . (3.36) 
which is about (3.3) times too big. The parabolic fit yields a well depth of 

D P3 (2.9 ± 0A)MeV (3.37) 

with the mass constraint 

H Q > (77 ± 6)GeV , (3.38) 

in order to obtain binding. 

Next the potential in equation 3.29 is extended to include all of the flux-bubble 
interactions, where the colour is fixed to each of the quarks: 

V = a ( r Q? ~ r o) #(r Qg - - r ) ± a s ^ A PiPj (— — J 6(r - r PiP] ) , (3.39) 

min{Qg} i<j VViPj r / 

with particle index pk e {Qi,cjj\i,j = 1,2}, such that k = 1,2,3,4, and SU C (2) 
colour factor 



A 



PiPi 



-f if (ft, Pi) C {(q h Q 3 )} 
\ if (PhPj) C {(qi,qj),(Qi,Qj)} 



(3.40) 
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Notice that V can be rewritten so that its extension to equation |3.29| is more appar- 
ent: 



V 



E 

min{Q<j} 



o- (r Q g - r ) 0{r Qq - r ) - - a, 



1_ 1_ 



#(r - r Q g) 



+ Q>s 



E 



PiPj emin{Qq} 



To, 



0{ro-r PiP .) , 



(3.41) 



where min{Qg} is the complement of the set min{Qg} , and therefore defining more 



precisely the sum in the exponent of the variational wave function, equation |3 . 1 2 . 
Figure [£7| shows the results of the Monte Carlo for this potential. 

The Monte Carlo results were checked against the analytic solution given by 
equation |3.30| at R = 5fm only, since the system no longer decouples at the origin. 
The results are summarized in table |3~~3 . 



Table 3.4: Monte Carlo (MC) vs. analytic results for a flux-bubble potential. 



Parameters 


Analytic 


MC @ R = 5fm 


E min (a,f3) (MeV) 


1527.07 


1527.66 ±0.05 


a 


1.74 


1.71 


P (fm- 1 ) 


1.37 


1.38 



An attempt was made to fit U(R) to the Morse potential |34 , |35 

U M (r) = U {e- 2/3(r - ro) - 2e- /3(r - ro) ) + U t 

with binding energy 

/3%c 



E 



-D + %u{{v + \) 



2J2n Q D 



and depth D = U , where u = (3c y 2D/hq . For this potential 

{(3hc) 2 



> 



8Un 



(3.42) 



(3.43) 



(3.44) 
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Q 2 Potential Obtained by Using exp(-(/Sr)°) with V qq (r) = flux-bubble. 
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R (fm) 

Figure 3.7: The Q 2 potential obtained using the flux-bubble interaction with 
pseudo-hydrogen wave function. 
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Table 3.5: Summary of parameters for the fits in figure 3.7 



Potential 


Parameters 


Pseudo- Morse 


(3 (/m- 1 ) 


r (fm) 


a 


1.34 ±0.20 


-0.63 ±0.19 


2.10 ±0.37 


or (MeV/fm) 


Uoo (MeV) 


r' (fm) 


40 ± 15 


1527.6 ±0.5 


0.088 ± 0.031 


Yukawa 


V (MeV) 


a (fm) 


Voo (MeV) 


9.16 ±0.42 


0.498 ±0.012 


1527.70 ±0.02 


Parabolic 


C (MeV/fm 2 ) 


r (fm) 


Vo (MeV) 


31.7 ± 1.4 


0.4042 ± 0.0044 


1524.3 ±0.1 


Morse 


P (fm- 1 ) 


r (fm) 


U (MeV) 


2.62 ±0.38 


0.333 ± 0.039 


3.09 ±0.73 


f/oo (MeV) 




1527.8 ±0.5 
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is the corresponding binding constraint. Table [T5] contains a summary of the fits 
for the various potentials. The Morse and parabolic fits give, 

(/iQWse > (U±±)GeV, (3.45) 

OgUaboiic £ (79±7)GeV (3.46) 

with well depths 

(D) Morse « (3.1±0.7)MeV, (3.47) 

( J D)pa rab onc « (2.8±0.1)MeV. (3.48) 

Although the Morse gives a lower value it is a rather dubious one since the fit did not 
describe the shape of the knee or the minimum of the potential (see Morse fitting 
insert in figure |3.7f) very well. Regardless, there is no binding. 

Referring to equation |3.20| and to table |3.5| yields an exchange mass of 

m cx « 0(396 ± 10)MeV , (3.49) 

which is (9(2.8) times too big. 

Finally the flux-bubble potential is extended to allow the colour to move around. 



The potential is similar to that of equation |3.39| except that now the particle indices, 
Pk, carry colour degrees of freedom, 

min{^g f } L 1 V <Mj '0/ 





3 




"I" 


i r 









'7j 



+ a s E a p iPj ( — r ( r o - r w ) > ( 3 - 5 °) 

PiPj emm{qiqj} 

where X PiPj is defined just below equation 3.1, and 

q k e {R i ,B i ,T U b i \i= 1,2} C SU C (2), (3.51) 



q k e {R i} B h f h b i \t = 1,2} C SU C (2), (3.52) 
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i.e., b ~ f,...,etc. & SU C {2) ~ SU C {2) , 

where the capital case letters are for the heavy quarks, the lower case letters are for 
the light quarks, the letters r and R represent the red quarks, and the letters b and 
B stand for the blue quarks. A further expansion of V leads to the more useful form 



V 



E 

min{qiqj} 



° (r mj - r ) 9(r q .g. - r ) - - a 



1 



1;, - 



3 

d s 

4 



E 



1 



qtqj e mm{qiqj} ^ 

E 



r 0/ 



8{r ~ r qigj ) , 



(3.53) 



where 



mmlqiqj} = mm{q i q j } + U min{g i g j }_ , 

plus, "+", means attractive, and minus, means repulsive. Figure 3J3 shows the 
results of the Monte Carlo for this potential. 

The Monte Carlo results were checked against the analytic solution given by 
equation |3.39| at R = 5fm oo). The results are summarized in table |3~^. 



Table 3.6: Monte Carlo (MC) vs. analytic results for a flux-bubble potential. 



Parameters 


Analytic 


MC @ R = 5fm 


E min {a,{3) (MeV) 


1527.07 


1527.62 ±0.04 


a 


1.74 


1.73 


P (fm- 1 ) 


1.37 


1.37 



Table 3./ gives a summary of the global and local fits to the potential. 



One of the most noticeable peculiarities of figure |3J| is the apparent linearity of 
the inside of the potential. In fact, a linear fit to 



y(r) = m(r — ro) + b 



(3.54) 
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Q 2 Potential Obtained by Using exp( — (/?r) a ) with V qq (r) = full — model. 
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Figure 3.8: The Qi potential obtained using the full flux-bubble interaction, in which the 
colour is allowed to move around, with a pseudo-hydrogen wave function. 
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Table 3.7: Summary of parameters for the fits in figure 3.8. 


Potential 


Parameters 


Linear 


m (MeV/fm) 


r (fm) 


b (MeV) 


910.0 ± 1.0 


-1.5706 ± 0.0047 


-534.6 ±4.2 


Pseudo- Morse 


(3 (/m- 1 ) 


r (fm) 


a 


0.18834 ±0.00033 


-4.423 ± 0.012 


11.29 ±0.093 


a (MeV/fm) 


U x (MeV) 


r' (fm) 


-1310.2 ±5.5 


1524.8 ±0.5 


0.6470 ±0.0011 


V 




0.9447 ±0.0024 



in this region (see the insert in figure [T8] and the fitted results in table |3.7|) yields a 
slope of m « 910.0MeV/ fm ! This corresponds to the linear potential U ~ aR. A 
more subtle feature is at the origin (see the "blow up" insert in figure |3.8|) where the 
potential starts to plummet to — oo. This region is due to the coulomb attraction 
between the two heavy quarks. Out near R = lfm there appears to be a barrier 
and beyond this no more structure. Therefore, as two mesons, each containing a 
heavy quark and a light anti-quark, are brought together from infinity they feel a 
repulsive force. When near enough, flux-tubes are exchanged, and the two mesons 
dissociate into one meson containing two heavy quarks and another containing two 
light anti-quarks. This situation shows that SU C (2) , with moving colour, does not 
make a good model of nuclear matter. 

Table [3i| gives a summary of all the properties of the Q2 potential for the linear, 
linear-plus-Coulomb, and flux-bubble (with fix colour) models that were studied in 
this section. In general the extensions of the basic linear potential model to 
include the colour- Coulomb interactions did not alter the Q 2 potential significantly 
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Table 3.8: Summary of the properties of U(R). 





Fig. 3.5 


Fig. 3.6 


Fig. |3 


D (MeV) 


3.4 ±0.1 


2.9 ±0.1 


2.8 ±0.1 


fi min (GeV) 


55 ±6 


77 ±6 


79 ±7 


m ex (MeV) 


638 ± 15 


459 ±4 


396 ± 10 



enough to achieve binding. In fact, it only exacerbated the problem. However, these 
extensions did lead to a slight softening of the potential which perhaps suggests 
that there might be more to meson exchange than just flux-tube swapping. When 
the colour was allowed to move around a rather unphysical situation occurred which 
suggested that there was perhaps a problem with using SU C (2) or with the variational 
wave function itself — perhaps even both. In the next section a more detailed study 
of the properties of the wave function is carried out. 



3.4.3 Back to the SU e {2) String-Flip Potential Model 

When the SUe(2) string-flip potential model was investigated || the parameter a in 



the variational wave function, given in equation |3.2j , was fixed by requiring that it 
minimize the total energy at zero density. Since this could be done analytically it 
allowed for a reduction in the number of variational parameters used in the Monte 
Carlo. It was assumed that the constraint would have very little affect on the physics 
as a function of density since the results varied by about 1% for 1.5 < a < 2.1 , 
at zero density. The validity of this claim can now be checked more thoroughly by 
using the Q2 mini-laboratory. 



Figure GO shows a plot of the Monte Carlo results, obtained via equations 3.10 



3.12| and |3.13| , for U(R) where a is allowed to vary, a = 2.00, and a = 1.74. 
a = 2.00 was the value used in the old SUe(2) model f|, and a ~ 1.74 was the value 
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Figure 3.9: U(R) where a is, allowed to vary, fixed at 2, and fixed at 1.74. These potentials 
were found by minimizing (U(R)) a n , for fixed values of R, using the distributed minimization 



algorithm (§ |3.4.l|) . The (U(R)) a ,p , were evaluated using the Metropolis algorithm (§ |2.4| ). 
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that minimized U (R) at zero separation. The values of U(R) at the end points of 
the curves, from R = Ofm and out to R = 5fm, were checked against the analytic 



solution given by equation 2.2S and have been tabulated in Table 3.9 



Table 3.9: Monte Carlo (MC) vs. analytic results for the SUg{2) ^-potential. 



Parameters 


Analytic 


MC @ R = Ofm 


MC @ R = 5/m 


E min (a,(3) (MeV) 


1709.61 


1709.42 ± 0.04 


1709.49 ±0.04 


a 


1.75 


1.74 


1.72 


P (fm- 1 ) 


1.37 


1.37 


1.38 


E min (2.00,(3) (MeV) 


1714.38 


1714.83 ±0.13 


1714.72 ±0.12 


(3 @ a = 2.00 (fm- 1 ) 


1.27 


1.25 


1.28 


£ wi „(1.74,/3) (MeV) 


1709.62 


1709.43 ± 0.04 


1709.43 ±0.04 


(3@a^ 1.74 (fm- 1 ) 


1.37 


1.38 


1.38 



Figure 3.10 shows the variations in a and (3 as a function of R. It would 
appear from this figure that when a is left free to vary, a and f3 become extremely 
correlated. Regardless, looking at the maximum fluctuations about the central values 
of the parameters a , (3 , and U(R) , from the information given in tables (j.lCj and 
[□J a w (1.88±0.13) , (3 w (1.32 ±0.05)/m, and (U(R)) R w (1710 ±A)MeV^ which 
is certainly less then a 1% effect in the energy 0. However for (3 it would appear to 
be an (9(4)% effect which may have a slight affect on the rate of nuclear swelling. 

Surprisingly the U(R)\ a= 2 gives a much deeper well than expected if a was left as 
a free parameter: i.e., D w 0(8)MeV , via table |3TT0| . However, this well is not deep 
enough to give binding: i.e., a parabolic approximation about the minimum requires 
> 0(15)GeV , via equation [3.26| and table |3.10| . For equal mass constituent 
quarks this is expected to be a much graver situation. 



1 i.e., (U(R)) R M (E mln (a = 2)\ R=00 - y )/2 ± ((U(R)) R - y )/2 where y is from table 3.1 
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Figure 3.10: Graph of a(R) & j3(R) for the SUe(2) Q 2 -potential model. The dashed lines 
are for a = 2.00 and a w 1.74. The noise is mainly due to shallow minima on the energy 
surfaces, U(R). The curve on the upper plot and the two upper curves on the lower plot have 
been parameterized by ae~ bH R c , which are provided to guide the eye; they have no physical 
significance. 
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Table 3.10: Summary of the parameters used in the fits for the a — 2 curve in figure B.9 
See equations 3.15| and 3.22 for the definitions of these parameters. 



Potential 


Parameters 


Parabolic 


C (MeV/fm 2 ) 


r (fm) 


Vo (MeV) 


49.9 ±3.5 


0.3944 ± 0.0046 


1706.7 ±0.1 


Pseudo-Morse 


ft (fm- 1 ) 


r (fm) 


a 


1.062 ±0.013 


-0.607 ±0.016 


2.768 ± 0.044 


a (MeV/fm) 


Uoo (MeV) 




39.68 ±0.87 


1714.70 ±0.04 



For certain fixed values of a , a ~ 1.74 in figure |3.9| for example, the potential 
gives a slight short range repulsion. For the range of a values being considered here 
it is well inside the potential. However, for values of a ~ 1.0 the effect becomes 
quite dramatic. 

In general, it is fairly safe to assume that the predicted outcome of the results 
of past papers [§, || will not change significantly if a is allowed to vary. However, 
because the wave function appears to be correlated in both a and ft , this seems to 
suggest that another wave function should be considered in these models. 



3.4.4 A New Wave Function 

Wave functions for the old SUg(2) string-flip potential models have been examined. 
For the analysis based on the Q2 system there appeared to be something patholog- 
ical about the wave function that was used in these models. Further, despite this 
shortcoming, it was thought that this would not have any significant effect on the 
outcome of the physics predicted by past models. However, the purpose of using 
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the Q2 system was not only to study the properties of this wave function but to 
ascertain the possibility of devising a wave function with very few parameters that 
would take into account the locality of the flux-bubble interactions. 

An interesting place to look for such a wave function is in the similarity between 
the Q2 mesonic-molecular system and H2 molecular system. Recall that the key 
reason for the atomic binding is the screening effect caused by the electrons which 
are for the most part "localized" in between the protons. This localization is achieved 
by using a variational wave function that was the superposition of the direct product 
of two ground state hydrogen atoms, |34], |38fl 



~P(r eiPl +r e2 p 2 ) -P(r e2Pl +r eiP2 ) 
^ ~ e + e . (3.55) 

Although the Q2 system is far removed from its H2 cousin from a dynamical point 
of view, and the motivations for achieving localization are quite different, this varia- 
tional wave function does solve the problem for the H2 system. Therefore, it would 
seem plausible to make the following ansatz: 

w a ,3 = e + e . (3.56) 

If qiQi and (I2Q2 represent two separate mesons then the first term represents the 
internal meson interactions while the second term represents the external meson 
interactions. Notice that the external interactions shut off as the separation, R , 
between the two heavy quarks becomes large, 

hm^(i?)=e -r(r?lQl+r?2Q2) (3.57) 

R— >oo 

which is the desired property. 

Using equation p.ll| and equation |3.57| the kinetic energy contribution to the 



effective potential U(R) in equation 3.13 now becomes, 

f = 2f_ s -F 2 (3.58) 
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— 23 a (r- 4- r- ) 

— V t L (« + 1) (rt Q \ + r-l) - a p (r-- 2 + r*£)] * ^ 



-2/3 Q (r| Ql +r| lQ2 ) 



+ [(« + 1) (^Q 2 + ^ 2 Q 2 + rgt + r££) 



— n>fl a (r 2a ~ 2 4- r 2a - 2 4- r 201 ' 2 4- r 2a ~ 2 ') 
a P VqiQl + T q 2 Q 2 +T Q2Ql +r qiQ 2 ) 



— R a (r a 4- r a 4- r a 4- r a \ \ 
xe P VqiQl ~+~ r q 2 Ql ~+~ r qiQ2 ~+~ ^Qj-jy^ \ 



(3.59) 



where 



-1 



£Vjln* 



(3.60) 



/(a + 1) 
4m \ * 



r.a-2 _|_ r a-2 
9lQl ~ <3 2 Q 2 



-^(^ 1 Q 1 +r| Q2 ) 



, , a-2 , o-2\ vmOi" 1 " qiQz 

+ Vg 2 Qi + r giQ 2 J e 



, [^2 , 2 _ r>2\ a-2 a-2 , / 2 , 2 _ e>2n a-2 a-5 

+ LVgiQi ~+~ r 5iQ 2 ^ J r giQi r qiQ 2 ~+~ V g 2 Qi ~+~ r ? 2 Q 2 n > r q 2 Qi r q 2 Q 



_oa( a _|_ a , a , a 
„2a-2 , „2a-2 , 2a-2 , 2a-2\l . 1 V 9lQl ^ qiQ\ ^ qiQ 2 ^ q 2 Q 2 



_ ( r 2a-2 , 2a-2 , 2a-2 , r 2a-2Ni 

V giQi + r ? 2 Q 2 + r q 2 Qi + r 5iQ 2 )\ 6 



(3.61) 
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and 



2m 



£(V ff ln¥) s 



2 o2a 



a 2 P 



2m q \ * 2 



C 2a-2 , r 2a-2N 



-2/3 a (r| Ql + r« lQ2 



(3.62) 



+ [( r 5iQi + r \xQi Rl ) r q x Qi r qiQi + ( r i 2 Qi + r ? 2 Q 2 ^) r ? 2 Qi r t 2 Q 



xe 



^ a ( r |iQi + r | 2 Qi + r 5iQ 2 + r 5 2 Q 2 



(3.63) 



It is interesting to note, and it also serves as a good check of the computation, that 
in the large R limit the kinetic energy corresponds exactly to the kinetic energy for 
the old wave function in the same limit: 

af3 a 



lim T(R) 

R^oo 



2m n 



g 2 y 2 



2m, 



5>a - wr, iQi r) + 1] 



(3.64) 



9 \i=l 



c/. equation |3.14| . Also direct evaluation of the RHS leads to (assuming \I/ is properly 
normalized) 

limf( J R) = ^/3 2 , (3.65) 



R—>oo 



which is just the kinetic term for the analytic solutions given in equations [3.16| and 

oq . 

The Monte Carlo computations that were done, using the pseudo-hydrogen wave 
function, (equation |3.12|) in § [3.4.2| , have been repeated here for the pseudo- 
hydrogen- molecular wave function, \& # 2 (equation |3.56| ), and are shown in figures |3~TT| 



CHAPTER 3. FLUX-BUBBLE MODELS AND MESONIC MOLECULES 



72 



through |3.14j. Table |3.3| contains a summary of the fits for these figures, table p. 12 



contains a summary of the checks done against the analytic results and table p. 13 
contains a summary of the properties of U(R) . 
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Q 2 Potential Obtained by Using an H 2 Wave Function with V qq (r) = linear. 
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Figure 3.11: The Q 2 potential obtained using a linear interaction with ^h 2 
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Figure 3.12: The Qi potential obtained using a linear-plus-Coulomb interaction with ^ 
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Q 2 Potential Obtained by Using an H 2 Wave Function, V qq (r) = flux-bubble. 



1540 - 



1 520 



500 



1480 - 



1 460 



440 - 



420 




Figure 3.13: The Q2 potential obtained using a flux-bubble interaction with *S?h 2 
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Q 2 Potential Obtained by Using an H 2 Wave Function, V qq (r) = full-mode 
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Figure 3.14: The Q2 potential obtained using the full flux-bubble interaction, in which the 
colour is allowed to move around, with ^ h 2 ■ It should be noted that there is a similar effect 



in the neighborhood of R = as shown in the "blow up" of figure 3.8 
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Table 3.11: Fits for figures pi] through [3l4| 



V 



Parameters ([Energy, Mass]~ MeV , [Length] ~ fm) 



Fig. BTT 



Fig.^TJ 



Fig.Q 



Fig. KU 



P 

s 
c 
u 

d 





M 



r 

s 
c 



0.727 ±0.060 



0.758 ±0.027 



0.736 ±0.056 



0.686 ±0.011 



r 



-0.957 ±0.096 



-0.893 ±0.057 



-0.94 ±0.12 



-1.527 ±0.022 



a 



2.50 ±0.11 



2.429 ± 0.062 



2.45 ±0.13 



1.367 ±0.034 



a 



308.3 ±8.9 



304.5 ± 8.9 



313 ± 11 



-260.8 ±9.3 



1.926 ±0.051 



1709.5 ±0.4 



1527.7 ±0.4 



1528.1 ±0.5 



1525.6 ±0.79 



Y 
u 

k 

a 
w 
a 



4400 ± 1700 



4500 ± 2000 



4000 ± 1800 



6000 ± 2200 



0.327 ±0.027 



0.335 ±0.031 



0.343 ± 0.032 



0.321 ±0.025 



V' 



1709.7 ±0.5 



1527.9 ±0.5 



1528.0 ±0.6 



1527.8 ±0.5 



P 

a 
r 

a, 

b 



C 



203 ± 11 



198 ± 11 



198 ± 10 



r 



0.591 ±0.047 



0.590 ± 0.048 



0.590 ± 0.048 



2/o 



1627.6 ±0.6 



1446.0 ±0.6 



1446.1 ±0.6 



1.536 ±0.013 



M 
o 
r 

s 
c 



r 



0.4613 ±0.0028 



Un 



85.47 ±0.69 



U n 



1535.0 ±0.6 



C 
o 
u 

1 



450.39 ±2.9 



2.628 ± 0.042 



1609.8 ±2.3 
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Table 3.12: MC vs. analytic results for figures [3. 11| through |3.l4 



Method 


R(fm) 


E min (MeV) 


a 


P (fm- 1 ) 


Analytic Eq. |3.16 


O/oo 


1709.61 


1.75 


1.37 


MC Fig. JJ.11 





1709.647 ±0.032 


1.77 


1.36 


5 


1709.485 ±0.033 


1.75 


1.36 


Analytic Eq. |3~30| 


O/oo 


1527.07 


1.74 


1.37 


MC Fig. p. 12 





1527.621 ±0.043 


1.73 


1.38 


5 


1527.651 ±0.057 


1.70 


1.38 


MC Fig. 3.13 


5 


1527.715 ±0.042 


1.73 


1.37 


MC Fig. 3. 14 


5 


1527.703 ±0.043 


1.73 


1.36 



Table 3.13: Summary of the properties of U(R) from parabola and Yukawa data in table 3.11 



with Uoc extracted from analytic results in table 3.10. For Morse data D w (85.47±0.69)Mey 



i^ mi „~ (134±2)M e y. 


Fig. 3.11 


Fig. 3.12 


Fig. 3.13 


Fig. 3.14 


D {MeV) 


Eq. p.24 


82.0 ±0.6 


81.1 ±0.6 


81.0 ±0.6 




/i min (MeV) 


Eq. p.26 


588 ± 33 


586 ± 34 


588 ±31 




m cx (MeV) 


Eq. p.20 


603 ± 50 


589 ± 55 


575 ± 54 


615 ±48 
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It can be immediately seen that there is a dramatic contrast between the figures 
for ^/h and ^h 2 ■ F° r figures |3.11| through |3.13| the wells are much deeper and 
the binding energy constraints have come down considerably, see table |3.13j for D 
and fi min (< hq). These potentials are now strong enough to bind the more massive 
quarks, such as c and b. This should be of no surprise as the adiabatic approximation 
used requires that the quarks be massive. The effective Yukawa masses, m cx in 
table |3.13| , are roughly the same order of magnitude as the case, table |3~8| , but 
with a slightly less dramatic softening effect. However, these masses are too large to 
explain pion exchange. 

An attempt was made to fit the flux-bubble model to a Morse potential but as 
before the fit failed (see insert in figure |3.13| ). It is interesting to note however that 
the pseudo-Morse potential of the form 



U(r) ~ {e-^-ro)" _2e-^-^ Q } ar 



(3.66) 



[cf. equation |3.15|) describes U(R) in figures |3.11| through |3.13| quite well. 

The final figure, figure |3.14| , of the flux-bubble model with moving colour is quite 



intriguing. The anomalies in figure |3.8| have disappeared; the light quarks have not 
drifted away as an isolated pair to leave a linear potential between the heavy quarks. 
In fact, there is no more linearity inside the well. Naively, this seems to suggest 
that the problem was with the wave function and not SU C (2). However this is not 
quite the case, since the old wave function gave an interior well depth twice as deep. 
Therefore, it is energetically more favourable for the Q2 system to dissociate into 
two isolated mesons; one with two light quarks and the other with two heavy ones. 

An attempt was also made to find the shape of the interior of the well, and it 
was found that it fitted to a Coulomb potential of the form, 

-07* 



V(r) = - 1 
r \ 



(3.67) 
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see insert in figure |37H| , and table pTTT] for fitted data), with strength a m — (150. 39± 



2.9)MeV fm . The term in the brackets was included to mimic the overlap between 
the charge distributions of two mesonic systems. In terms of a s , the hyperfine 
constant for this region of the potential is 

a Q2 « (7.62±0.15)a s . (3.68) 

The interior part of the potential is quite deep, O(270)MeV , and bottoms out at 
(9(1255) MeV, at which point the —a s /r term for the heavy quarks kicks in (i.e., 
at R = O.lfm). The exterior part of the potential fits to a Yukawa potential with 
m ox « O{600)MeV, via table 

3.5 Discussion 

Various aspects of model building for nuclear matter have been examined in the Q2 
system. The ramifications of these investigations in regards to nuclear modeling will 
now be discussed. 

Perhaps the most important result is the new variational wave function, equa- 
tion |3.56| . This new wave function made a large change in the depth of the Q2 
potential well. The depth increased by a factor of 27, deep enough to bind heavy 
quarks: i.e., m q > 0(m c ). The wave function also fulfills the requirement of handling 
local flux-bubble interactions, which becomes apparent when looking at the before 
and after pictures (of moving colour) in figures |3.8| and P . 14] , respectively. Therefore, 
in SUg(2) for a many quark system this would suggest the following ansatz: 

*~Perm|* H (r PiP ,)| II M r P J\ ( 3 - 69 ) 

colour 

where 

**(r w ) = e ^ , (3.70) 
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Perm|^ij(r PiPj .)| is a totally symmetric pseudo- hydrogen wave function and |$(r Pfe )| is 
a totally antisymmetric Slater wave function (for the full three quark system a similar 
wave function would apply). This does not necessarily mean that this wave function 
would lead to nuclear binding, since binding was only achieved for relatively very 
heavy quarks in the Q 2 system. However, given the order of magnitude of increase 
in the Q 2 well depth it would seem quite plausible that it might be a strong enough 
effect to produce a shallow well in the nuclear-binding-energy curve. A simple test 
would be to consider SUe(2) with just a string-flip potential, with a fixed, in which 
case scaling is restored and the Monte Carlo becomes quite straightforward to do. 

The flux-bubble model proved quite successful at combining the colour- Coulomb 
interactions with the flux-tube interactions. Although these interactions, in general, 
had very little affect on the Q 2 system it was useful in demonstrating that extending 
the flux-tube model to include local perturbative interactions can be done. Further- 
more, it was not surprising that this had very little effect as it was hypothesized, in 
chapter 2, that the hyperfine interactions should play the dominant role. Therefore, 
it would prove most interesting to investigate the effect of adding more perturbative 
interactions to the Q 2 system. With the addition of SU C (3) this would lead to a more 
realistic model of mesonic molecules which perhaps could be tested in the laboratory. 

When the SU C {2) flux-bubble model was considered with moving colour the re- 
sults were quite interesting; the Q 2 system dissociated into one light and one heavy 
meson. However, this system is not physical. Perhaps a more useful model would be 
to consider the heavy quarks as a composite of two light quarks and using S77 c (3) 
instead. In this case it would not be possible for a flux-tube to form between the 
two heavy quarks causing the system to dissociate. 

Curiously if the heavy quarks were considered to be a composite of two light u 
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quarks then 

A* Qa{w} /M-»«O(0.6) (3.71) 

which is getting closer to binding, but is not quite sufficient. This would seem to 
suggest that perhaps slight nucleon deformations would help to enhance binding. 
Therefore, an interesting possibility would be to consider a many-body SUi(2) or 
SU C (3) (with moving colour) model in which there is an imbalance between the 
quark and anti-quark masses. In the case of SU C (3) , a strong enhancement would 
be expected given the depths of the wells observed in the Q2 system with SU C {2) . 

The distributed minimization algorithm proved to be quite successful. The algo- 
rithm as it stands should have very little difficulty handling SUe(2) models in which 
both a and f3 are allowed to vary. The reason this can be said with any confidence is 
that the energy- surface topologies of the many-body SUe(2) models are very similar 
to those in the Q2 system. From a practical point of view the projected time for a 
network of m computers is 

T~0(—)hrs, (3.72) 
m 

for a fixed density parameter, p, which is just barely tolerable for m = 8. If a is 
fixed then 

162 

r~C>( )hrs, (3.73) 

m 

which is more reasonable. 

For SUe(3) the projected times increase by a factor of 0(4). The simplest way to 
circumvent this problem is to increase the number of computers. With the resources 
here, at Carleton, it is possible to go to m = 16 by running off of the 8-node CPU 
farm with some slight modifications to the client-server routine that distributes the 
jobs: i.e., they have to be modified to use "unused" machine time. 

Neural networks were also investigated as another possibility of decreasing CPU 
time. Some preliminary work was done using a backpropagation network (BPN) |5{J 
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to perform flux-tube sorting. Within about 5 minutes of training, the BPN found the 
minimum of the string-flip potential to within 20%. The BPN was then compared 
to the fragmentation procedure (for a 7-quark qq system) for several hundred flux- 
tube configurations and it obtained a better solution about 1% of the time. The 
BPN method sorts in constant time and therefore is not dependent on the density, 
p, and has the potential of handling a much larger number of quarks. However, at 
its current stage of development the projected learning rate required to obtain an 
accuracy of 1% is about two months. Nonetheless, the work done with BPN's thus 
far seems to indicate that they hold a very promising future. 

The program for finding a model of nuclear matter is now quite clear: examine 
the effect of 

• the new wave function on the old SUe(2) model ||, 

• SU C (3) on the Q 2 system, 

• adding perturbative interactions to the flux-bubble model on the Q 2 system, 

• the interplay between the aforementioned scenarios, 

• the most relevant of these scenarios on the many-body SUe(2) model, 

• m q rrig on the many-body SUt(2) model, 
and finally 

• move on to studying SUi(3). 

All of these investigations, with perhaps the exception of the last, are now possible 
to perform with the tools that were developed in this thesis. 
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3.6 Conclusions 

The Qi system has proven to be a very useful aid for trying to sort out the complex- 
ities of model building for nuclear matter. The details of the mechanics, from wave 
functions to dynamics to practical computing methods, of the flux-bubble model have 
now been thoroughly investigated. It appears that the flux-bubble model may prove 
to be very successful, not only for modeling nuclear matter but also for modeling 
mesonic molecules as well. 



Chapter 4 



L+L Production in E 



6 



4.1 Introduction 

In this chapter the production of heavy leptons pairs, L + L~ , via superstring inspired 



E 6 models, at high energy hadron colliders will be investigated [40]. Here several 
underlying assumptions will be made about E 6 models in order to restrict the com- 
putation to a manageable, but reasonable, calculation. The specific model will be 
presented followed by a summary of the calculation and then its phenomenological 
consequences. 

NOTE: Many aspects of the rank-5 models, that will be considered here, are cov- 
ered in the literature. Unfortunately, when trying to extract the particular model 
dependent information needed for L + L~ production, it appeared that the existing 
literature was not consistent. Therefore, it was felt that in order to avoid any am- 
biguities, the model should be carefully reconstructed from the ground up. When 
constructing the model, careful attention was paid to being as consistent as possible 
with the literature concerning; factors of two, hypercharge conventions, signs, am- 
biguous notational subtleties, etc. Much of the analysis of the model was done by 
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using Mathematica jj^ffl to generate the various couplings, mass matrices, etc., di- 
rectly from the superpotential. This enabled easy comparison with various literature 
sources ]73| , \Th\ [Q , \4% , \jphJ . Differing conventions and normalizations aside, the 
most significant problem arose with the charged- Higgs, equation 4-3% and pseudo- 



scalar-Higgs, equation \4 . 33j , mass terms; a factor of two was missing in front of the 
sin p cos (3 terms, op. cit. For example, in the case of the pseudo-scalar-Higgs the 
aforementioned authors disagree by an overall factor of two in their mass-mixing ma- 
trices but not in their eigenvalues. As a result, the analysis of the mass constraints 
in the Higgs sector J^/ had to be re-evaluated, figures \4-3j through \j.b\ . In addition, 
Appendix [£J contains a summary of the couplings used for L + L~ production which, 
in general, could not be obtained from the literature. 

4.2 A Low Energy E 6 Model 

In chapter 1 a general overview of superstring inspired E 6 models was given. In 
addition, several comments were made about the type of model that would be pre- 
sented. We shall now expand on these assumptions to determine their low energy 
consequences. 

There are many ways of breaking E 6 down to SM energies. Invariably these 
breaking schemes lead to SM phenomenologies which contain extra gauge bosons. 
Here a rather simple model was chosen in which only an extra Z, the Z' , is produced 
|,§: 

E 6 — > SU(3) c ®SU(2) l ®U(1)y®U(1) Ye 



(c/. figure |1.6Q . In general the Z' can mix with the SM Z to produce the mixed 
states 

Z\ = cos Z + sin Z' , (4-1) 



CHAPTER 4. L+L- PRODUCTION IN E, 



87 



sin Z + cos Z' , 



(4.2) 



(cf. equation |4.45| .) 

Recall that in order to avoid potential problems with flavour-changing-neutral 
currents, at the tree level, a basis was chosen in which the third generation of primed- 



exotic-sleptons were assigned to play the role of the Higgses fll3| , |15 




(4.3) 



(i.e., R3 = L 3 ) or by redefining L 3 , R' 3 and v 1 !^ , in terms of the complex-isodoublet 
fields, $1 and <3> 2 , and the complex-isoscalar field, $ 3 , respectively, equation |Q 
becomes 



$1 



This assignment was accomplished by setting 



4>V 




( A+\ 





















(4.4) 



\333 _ \3j 
^4 — ^4 



3i3 
k 



= \p = Af 1 = 



\{ 3h 



Af 3 ± 




(4.5) 



in the superpotential, equation |1.2| , where the ijk 1 s are generation indices. Therefore 
in order to avoid lepton-number violation the lepton-numbers of all of the primed 
and double-primed exotic-leptons must be zero. 

Further restrictions were placed on the superpotential by requiring that the 
baryon and lepton numbers of the exotic-quarks, q' , of the 27's (figure |i.5p, are 



the same as those of their non-exotic SM counterparts Jj]| : 

\ = A7 = As = Ag = A10 = . 
Also An was set equal to zero in order to avoid any fine tuning problems with the 



B(q' L ) = i L(<Q = 



vf masses 13 
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With all the aforementioned assumptions about the Yukawa couplings, the su- 
perpotential now simplifies to 



(4.6) 



where the A's were chosen to be real, plus similar terms for the other generations 
and their cross-terms. The $^ = ^{ipAiA) are the superfields which contain a 
two-component-spinor field, ipA , and a complex-scalar-singlet field, A . Table [4.1| 
summarizes the particle properties of this model. 



Table 4.1: Table of SU(3) C ® SU(2) L ® U(1) Y ® U(1) Ye particle properties. 



4 C 



L' 



R' 



1/2 

-1/2 



1/2 
-1/2 



1/2 
-1/2 



1/2 
-1/2 



Y 



1/3 



-1 



-4/3 



2/3 



-2/3 



2/3 



-1 



2/3 



-1/3 



2/3 



-1/3 



2/3 



5/3 



-4/3 



-1/3 



-1/3 



-4/3 



5/3 



Q 



2/3 
-1/3 



-2/3 



1/3 



-1/3 



1/3 



D 



1/3 



-1/3 



■1/3 



1/3 



-1/3 



-1 



-1 



The superpotential specifies all of the couplings between the particles of the 27's. 
According to appendix B of Haber and Kane [[EJ , the Yukawa interactions are given 
by 



Yuk 



1 

2 



' d 2 W 
dAidA, 



i>'s = 



' d 2 W 



■<P's = 



(4.7) 



and the scalar interactions are given by 



V = V F + V D + Vs of t 



(4.8) 
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In equation p~8|, 



Vf 
V D 



F*F, 



% 1 1 



[D a D a + (£>') 



/\2l 



(4.9) 
(4.10) 



with 



D' 



gA*T^Aj . 
\g%A*A 



■4>'s = o 



(4.11) 

(4.12) 
(4.13) 



g represents an SU(N) coupling constant with generators T a , and g' represents a U(l) 
coupling constant with hypercharge Y. Vsoft is a soft-SUSY-breaking term that was 
put in by hand in order to lift the supersymmetric-mass-degeneracy between m^ A 
and rriA '■ 

Vsoft = V, + V I+ V Hsoft} (4.14) 

with 



2 Re [\ 1 A u iQ T T 2 $ 2 u$ L + X 2 A d i^r 2 Qd* R + X 5 A df d' R d' L ^ 3 } , (4.15) 
Vj D M 2 L \L\ 2 + M e 2 e* R e R + M^ e * R £ eR + 2X 3 A e Re [*$fr 2 Lg R ] , (4.16) 
V Hsoft = ^|$ 1 | 2 + ^|$ 2 | 2 + ^|^3r-^AA(i$[r 2 $2^3 + /i.c). (4-17) 



The coefficients M\ , Aa , /i 2 , and A are the soft-SUSY-breaking parameters, and 
A = Af 3 . 

The Higgs potential can be extracted directly from equation and is given by 



13, 43, 44 



V H = ^ So/t + A 2 (|$ 1 | 2 |$ 2 | 2 + |$ 1 | 2 |$ 3 | 2 + |$ 2 | 2 |$ 3 r) 
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+W + /xi$ii 2 - i^ 2 i 2 ) 2 + y\m 2 + 4i$ 2 r - 5|$ 3 | 2 ) 2 

+ {\g 2 - \ 2 )\^ 2 \\ (4.18) 



where g , g' , and g" are the SU(2)l , U(l)y , and U(l)y B coupling constants, respec- 
tively. The minimization condition JTB] 



0. 



(4.19) 



VEV's 



can be used to fix the /i 2 terms in Vn Soft [Eq. (|4.17|) 1, where the vacuum expectation 
values, V^-EVs, are given by 

' " AO 



Vi 

-4= if* 



otherwise 



(4.20) 



An analytical solution for the fi^s can be found; using Mathematica [B6 



A*? 
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/> 2 + - ^We., + ^(g 2 + g' M - 4Wi 



£ ijk 

j<k 



{Vj+Vk) A 2_ VjU k 



2 Avi 

where = \£ijk\ ■ The kinetic terms for the scalar fields are given by [|6 



Ck.e. 2 



2 = l(^-|^)^| 2 , 



(4.21) 



(4.22) 



with (cf. [T 

Gfi = (gr 3 sin 8 W + g'Y cos 6^)^ + (# r 3 cos - ^ sin 9 w )Z,j, 

+V2^[r + tr;- + r^W+] + , (4.23) 

where t± = ~(ti ± «r 2 ), Tj|$ 3 ) = 0, and </ ~ |13], 0. The Tj's are the Pauli 
matrices acting in isospin space. 

The <&j fields have complex components, 0" , which were chosen to be of the form 



m m (of. 0) 



i 

71 



(4.24) 
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where <p a iR /\f2 and <f>yV2 are the real and imaginary parts, respectively. Therefore, 
the fields have a total of 10 degrees of freedom: four are "eaten" to give masses 
to the W , Z, and Z' bosons, and the remainder yield two charged-Higgs bosons, 
H ± , one pseudo-scalar-Higgs boson, P° , and three scalar-Higgs bosons, H® =1 2 3 . The 
mass terms for the Higgs fields can be obtained from the second order terms of the 



expansion of V#(0fc) about its minimum [fHfl , 
where 



d 2 V, 



H 



(4.25) 



(4.26) 



VEV's 



is the Higgs-mass-mixing matrix. Therefore the mass terms for the Higgs fields are 
simply 



-M 



D 



M 2 



P o 




31/ 



1 

2 



^IR - "l <P°2R - Vl 4>l R 



v 3 )M 2 h9 



( A.0 

°1R 
^R 





3R 



"3 J 



(4.27) 



The mass-mixing matrices are 



M 2 



V2 

[\g 2 - A> 12 + \Au 3 {\g 2 



\ 2 )v 2 + \A — 



V_23 



(4.28) 



M 2 



\Av 3 



P o 



( 


V\ 


1 


V\ 






V2 


V-i 






1 




"2 








V\ 


"3 






V\ 


"2 


V\2 




\ 


V3 




v 3 


) 



(4.29) 
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M 2 



H9 



B x v\ + XA — B 2 v l2 - XAv 3 B 3 v n - XAu 2 
B 2 u 21 - XAu 3 B±v\ + \A — B 5 u 23 - XAi/x 
-63^31 — XAu 2 B 5 is 32 — XAui B§v\ + XA — 



where Uij = ViVj . In equation |4.30 



1 '" 2 + /) + ha" 2 B 2 = 2A 2 + y> 2 -\{g 2 + /) S 3 = 2A 2 - ^ 



5 ji 2 



B A 



r(<? 2 + </) + §<?" 



5 5 = 2A 



2 10 „// 2 

9 y 



r _ 25 „// 2 



(4.30) 



(4.31) 



The physical states are obtained by diagonalizing the terms in Cm . The eigenvectors 
for the charged and pseudo-scalar Higgs terms are respectively, 



G ± 



cos (3 (j) 2 + sin (3 (f)f , 
sin (3 4>f + cos /3 0^ , 



(4.32) 
(4.33) 



and 



g° 



G° 



where 



AAz/ 3 



\ 2m 



2 




2 ™3/ 



z/ 2 z/ 3 



(^ + ^)K 2 + ^I) L^2 



^3 ,0 Vl 



J 21 



1 + 



^3 



(4.34) 
(4.35) 
(4.36) 



v 2 = v\ + 1/| 1 an d tan /? = v 2 jv\. Here, the if are the charged- 



Higgs states with masses 



m 2 H± 



XAi>z 



sin(2/5) 

the P° is the pseudo-scalar-Higgs state with mass 



1 - 2 ^2 J m w > 



mpo 



AA^ 3 



1 + 



sin(2/3) r ' 4z/| 



sin 2 (2/3) , 



(4.37) 



(4.38) 
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and the G and G® 2 are the Goldstone-boson states with zero mass. The scalar- 



Higgs term can be diagonalized exactly using Mathematica p6fl , however the result is, 
in general, not very enlightening. For our purposes it suffices to resort to numerical 
techniques. In the unitary gauge (U-gauge) the Goldstone modes vanish, i.e., the 
G's = , and the fields become physical. This allows the change of basis: 



<f>t 



sin PH* , 
^23 P°, 

KV 13 P°, 
3=1 



(4.39) 
(4.40) 
(4.41) 
(4.42) 
(4.43) 
(4.44) 



where k = 1/ y v\v\ + v 2 v | . The Uij are the elements of the inverse of the matrix 
that was used in the similarity transformation to diagonalized the scalar- Higgs-mass 
term. With these transformations at hand it is now a straightforward matter to 
obtain all of the masses and couplings for the various particles in this model. 

The mass terms for the gauge fields can be found by transforming the kinetic 



terms for the <3>i fields, equation [4.22| , to the U-gauge basis, equations |4.39| through 
0|, yielding: 



C2e. ^ rn 2 w W+W- 



As a consequence, the W mass is 



-{Z Z'),M 2 Z _ Z , 



\ Z 'J 



(4.45) 



m 



w 



■ 2 2 

;9 " , 



(4.46) 
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and the Z — Z'-mass-mixing matrix is |13], [IT], [12], [|3 

m% 5m 



M 2 Z _ Z , 



Sm 2 m\, 



(4.47) 



with matrix elements: 



rrir 



m 2 z , 



5m 



^ 2 + </V, 



12 



g 2 + g"g"(Au 2 -u( 



(4.48) 
(4.49) 
(4.50) 



Diagonalization of M. 2 Z _ Z , yields the mass eigenstates given by equations [O] and [12 
with eigenvalues 



m 



Zi 



m 



m 2 z cos 2 (p + 5m 2 sin (20) + m 2 z , sin 2 
m 2 z sin 2 — 5m 2 sin(20) + m|, cos 2 



and mixing angle 



tan(20) 



26m 2 



m 2 z — m 2 Z i 



(4.51) 
(4.52) 

(4.53) 



Notice that in the large z/3 limit — > tt/2 and therefore Z\ — > Z and Z2 — > -Z 7 . 
In fact, for the range of VEVs that will be considered here m Zl < m Zi . Therefore, 
Z\ will be designated the role of the observed Z at facilities such as LEP or SLC. 
The mass terms for the fermions, and hence the Yukawa couplings, can be found 



by evaluating H Yu k 5 equation P~T| , in the U-gauge basis and then using Appendix A 
of Haber and Kane [45], to convert to four component spinor notation.0 The result 



is 



£y u k — — {^ lh>2 uu ^ 2l>1 + ^s^ 1 ee + -^ 4Z/ 3 e ' e> + -^5^3 dfd'} ■ (4.54) 



1 cf. equation C.28. For a more explicit example see section 4.2 of Gunion and Haber 
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Therefore, the Yukawa couplings for the first generation are given by: 



Ai 
A 2 
A 3 



m w sin (3 

gm d 

m w cos (5 
gm e 

m w cos (5 

V2 



A, 



"a 

"3 



■m d i 



(4.55) 
(4.56) 
(4.57) 

(4.58) 
(4.59) 



and similarly for the other generations. 

The sfermion masses are obtained by evaluating the scalar-interaction potential, 
V (equation fP|), and then transforming it to the U-gauge basis: 



C M 2 -V D -(fir R )M 




with 



M) 



M 2 LL Ml R 
M 2 LR M 2 RR 



(4.60) 



(4.61) 



being the sfermion-mass-mixing matrix. The mass-mixing matrix elements are found 



to be, using Mathematica |p6| : 

Mft = M2 + m^ + i(3-4x w )m|cos(2/3)-^'' 2 K 2 + 4z, 2 2 -5z, 3 2 ),(4.62) 



M, 



(«) 2 



M, 



RR 



LR 



M 2 U +m 2 u + -x w m 2 z cos(2/3) - g" 2 {v\ + 4z/ 2 - .V 
m u (A u - m e > cot/?) , 



(4.63) 
(4.64) 



for the UL t R squarks; 



M L f = Ml + m 2 - 1 (3 - 2x w ) m% cos(2/3) - 1 g" 2 [v\ + 4z, 2 - 5^ 3 2 ) , (4.65) 



M 



M, 



RR 

{df 
LR 



+m 2 d --x w m i z cos(2/3) + — (z/ 2 + Au 2 - 5u 2 3j , 
m d (A d - m e > tan/3) , 



(4.66) 
(4.67) 
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for the di,R squarks; 



~ 2 X X 

M% )2 = M\ + m\, + - m| cos(2/3) + — </' 2 (z/ 2 + 4z/ 2 - 5u 2 3 ) , 



M 



M 



id') 2 
RR 

(d'f 
LR 



1 



Mi +mi,-- x w m 2 z cos(2/3) + — g" {v{ + Au 2 - 5z/ 2 ) , 



72 



^12- 



(4.68) 
(4.69) 
(4.70) 



for the d' LR squarks; 



M 



RR 



M 2 e +m 2 e -x w m| cos(2/3) - ^ <?" 2 [v\ + 4z/ 2 - 5z/ 2 ) , 



m lr = m e (A e - m e > tan/3) , 



(4.72) 
(4.73) 



for the em sleptons; 



M 



LL 



Ml + ml + - m| cos(2/3) + — g" 2 (y\ + Av\ - 5i/ ; 



72 



M 



5 



M 



LR 



Ml-^g" iy\ + Av\ — 5z/ 2 ) , 
0, 



'37 > 



(4.74) 

(4.75) 
(4.76) 



for the z/ i jR sleptons, and similarly for the other generations. The mass eigenstates 
are given by 



/ A 



t 



cos Of sin 9f 
-sin9j cos9j 



( f 

JL 



f. 



R 



with mass eigenvalues 



(4.77) 



m 



h 



m h = 



M\ L cos 2 Of + M\ R sm(2^-) + M 2 RR sin 2 6; , 
M 2 LL sin 2 6 f - - M 2 LR sm(29 f ) + M 2 RR cos 2 6 f , 



(4.78) 
(4.79) 



and mixing angle 



M 2 
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Notice that for fairly large v% 

tan(20 / ) ~ O (^p) > ( 4 - 81 ) 

where the soft terms have been assumed to be large and degenerate. Therefore, in 
general, the mixing is only expected to effect the sfermions that have fairly heavy 
fermion partners. 

The supersymmetric partner, or spartner, degrees of freedom for the neutral- 
Higgs fields (neutral-Higgsinos) , 1/ , 1/ and i/£ , along with the spartner degrees 
of freedom for the neutral gauge fields (neutral-gauginos) , 7 , Z , and Z' , mix to form 
a (6x6) neutralino, x° > mass-mixing matrix.! Similarly the charged-Higgsinos, t' l 
and t'I , and the charged-gauginos, W , form a (2x2) chargino, x , mass-mixing 
matrix.! By virtue of supersymmetry the neutralino and chargino mass-mixing ma- 
trices contain the same Yukawa and gauge couplings as their spartners, modulo soft 
terms. 

The real and imaginary parts of the sfermion fields, i>' eL , v*° L , v"° L , v' h , u£ L , 
and u'' c L , yield two separate (6x6)-mass-mixing matrices for the neutral- unHiggses 
(cf. equation |4.27| for cf) iR and (j)^), which contain the Yukawa couplings given in 



equation [4.5| . In general, these mass-mixing matrices are expected to lead to very 



massive unHiggs states 15 



The spartner degrees of freedom for the neutral-unHiggses, v' eL , v'l L , v'£ L , v'^ , 
iC L , and u" c L , form a (6x6)-mass-mixing matrix for the neutral-unHiggsinos. There- 
fore, the neutral-unHiggsino mass-mixing matrix contains the same Yukawa cou- 
plings as their neutral-unHiggs partners. 

The sfermion fields e' L , e£ , jX' L , and , yield two separate (2x2)-mass-mixing 



2 A detailed study of the x° mass spectrum can be found in p8| . 

3 The full form of these mass matrices can be found in Ellis, et al., E3] and the details of how 
to obtain them can be found in appendix B of Haber and Kane . 
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matrices for the charged- unHiggses (cf. equation [4.27| for <pf). These matrices have 
a large number of unknown parameters and quite naturally acquire a very large mass 
(cf- tH). 

Finally, the spartner degrees of freedom for the charged-unHiggses give diago- 
nalized mass eigenstates (see equation [4.54|) which correspond to the charged heavy 
leptons. 



4.3 L + L Production Cross-Section 



Figure PO shows the Feynman diagrams used for computing the parton level gluon- 



gluon fusion to heavy leptons matrix elements. It was shown |49j that the E 6 matrix 




b) 




H»,P° 




Figure 4.1: Feynman diagrams for gluon-gluon fusion to heavy charged leptons. 

element computations are very similar to the corresponding MSSM calculation by 
Cieza Montalvo, et al., [0 and can be easily extracted from their paper. The matrix 
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elements are as follows: 



1. For the Z 12 exchange diagram shown in figure |4.1| .a 

Pl\ hc L L ±Z '-C L R ±Zl ] Z Z *(s) EiC q L Zi -C 9 R Zi )(l + 2X q l 



a L± 



a 2 a 2 s m\ 1 
YI%txx 2 w m^r , i 



where the left-right Fermion, / , couplings are given by 



U L,R 

r<f z 2 



cos <p sm i 
- sin 6 cos < 



/ 



r fz> 

\ L,R 



qj i 
(4.82) 



(4.83) 



such that 



r fz 
r fz' 



Ts L<R e f x W 
2 U 



(4.84) 
(4.85) 



where T 3r = —T£ L is the isospin, y'j R = —y'fc is the Y^-hyper charge, and e/ is 
the electric charge. 



2. For the Hf 23 and P° exchange diagrams shown in figure [4.1| .b 



(7 



9 P° 



a 2 a 2 s m 2 L 



h{K L±p0 ) 2 \C P \s)\ 2 \Y.K« p °m 2 q l 



2 

91 ' 



(4.86) 
(4.87) 



where the couplings K-f H i and K^ p ° are given by equations |C.29| through |C3S . 



3. For the H® 2 3 exchange diagrams shown in figures |4.1| .c and [4.1| .d 



- qHf's 



a 2 a 2 m L p L 



oiznx w {i x w ) - =1 - fc=1 



2 

<?fe 7 1 ' 



(4.88) 
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where the sfermion mass eigenstates, /i^ , couplings to the Hf are given by 

k[ H ^ = K f L f' cos 2 6 f - + K[f i sm29f + sm 2 6 fj (4.89) 
k[ H ^ = sin 2 Of-K^ sin 29; + cos 2 Of, (4.90) 

fH° 

where K A g (A, B = L,R) are the corresponding couplings for the sfermion 

fH° 

helicity states, jl,r, with mixing angle 9j. The K AB ' couplings are given by 
equations |(J.39| through |C.54 . 



4. For the q(Zi j2 -P°) interference terms, via figures LI. a and 4.1.b 

2 



a 



= ^^^^ ^L ± P0 p o z ^ ±Zi _^ ±Zi 

128?^ r</ v / L l a 



i=l 



X 



E E - <?£*] (1 + 2 ViJ) } • (4-91) 



5. For the (q — q)H® 23 interference terms, via figures O.b through 4~T.d 



a 



(q-q)Hf's 
L± 



-a 2 a 2 



2567nr 2 y (l - x w ) 2 



x 



E^°^[2 + (4\ - 1)/ 9 ] E^ L± ^c^(s) 

9 j=l 



*EE^f% + 2 V|)}- 

g fc=l 



(4.92) 



In the aforementioned list of cross-section equations: 

(1 — x)x 



1 dx 
— In 

o x 



1 



A, 



sin 



i I i 



2,/A r 



if A p > \ 



\ In 2 (S) - £ + «r In (^) if A p < I 



(4.93) 



is the loop function [16], where 
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and r± = 1 ± Jl — 4A P such that pe {f, /} 



c 



rrtr 



s — m|. + i m z .T z . ' 



m |fO pO ± i m }{0 pO-TpO pO 



(4.95) 
(4.96) 



are related to the vector and scalar propagators, respectively, where the widths, i.e., 
the IY^'s, are summarized in § |C.2| , and 



1- 



4m| 



(4.97) 



is the lepton velocity in the cm. frame. The details of the various components that 
have gone into this computation can be found in appendix y. Before the parton 
level cross-section can be used to compute the heavy lepton production rates some 
assumptions about the parameters and masses in the model must be made. 

The first thing that has to be constrained are the VEV 's. It is reasonable 
to assume that t>i/i>2^1, since m& << m t , for any reasonable range of Yukawa 



couplings [42, 44|. Now the ratios t>i/t>2 and V3/V2 can be constrained by looking at 
how the variation in the Z\ (i.e., the U Z") mass affects x w 



sin 2 6 W ) such that 



42 



2 

■ 2 a — 1 m W 
sin 9 W = 1 — 



< sin 2 6 



9' 



m 



w 



9 2 + 9" 



(4.98) 



where x w and calculated using the on-shell and MS schemes plf, respec- 



tively, ^^(m^) can be found by evolving ^^(m^) m 0.2319 ± 0.0005 [^TJ down to 



m 



p| , which gives i^fm^) ~ 0.233 ± 0.035 , with a 1 \ 



rrir 



127.9 ±0.1 51 



m 



z « (91.187±0.007)GeV |5T[, and m w « (80.23 ± 0.18)GeV J5J. Therefore given 



x w « 0.2247 ±0.0019 



yields 



A 



.7; 



.7; 



ir 



0.008 ±0.035. 



(4.99) 
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Figure 4.2: Plot of m Za and A contour lines as a function of V3/V2 and V1/V2 ■ The A contour 
lines are shown at the Oct , la , and 2a levels. The arrows point toward the allowed regions 
on the plot (c/. p2[ ). The m Z2 = 500GeV line shows the CDF and D0 constraints, assuming 
standard couplings [M. 
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Figure [4.2| shows the A contour line as a function of v 3 / v 2 and v 1 jv 2 , along with 
its la and 2cr level contour lines. Also shown are the m Z2 contour lines. Taking 
a la level constraint implies v 3 /v 2 > 0(3.5) (c/. @, 0) and m Z2 >O(200)GeV. 
Unfortunately these constraints are not that tight due to the large uncertainty in 
a(m z ). A stronger constraint can be found by using the CDF and D$ limits on 
the m Z2 mass [ 50 1 , assuming SM-like couplings, figure |4.2j . This constraint is fairly 
reasonable since Y E 's ~ C?(Y)'s (cf. table |4.1|). With these constraints v 3 /v 2 £ (9(7.5) 
and m Z2 > O(500)GeV. 

Figures |4l| through ^75] show mass contour plots as a function of m p0 and 
m H± for (v 1 /v 2 ,v 3 /v 2 ) = (0.02,6.7), (v 1 /v 2 ,v 3 /v 2 ) = (0.5,7.7), and (v 1 /v 2 ,v 3 /v 2 ) = 
(0.9,9.1), respectively, such that m Z2 lies roughly around the CDF and Dtf) limits. 
These figures are a fairly good representation of the behavior of the m„ contour 
lines as a function of V\jv 2 . For fixed V\jvi the contour lines change very little (i.e., 
<0(5)%) for O(10) >v 3 /v 2 ^j 0(4.5) . This region corresponds to the (vi/v2,v 3 /v 2 ) 
parameter space for mz 2 £ O(?>00)GeV depicted in figure Further examination of 
the other scalar-Higgses shows m H o is fairly insensitive to variations in m p0 and m H ± 
for fixed m Z2 and is slightly sensitive to variations in v 3 /v 2 , whereas the behaviour 
of ttItto appears to be quite sensitive to any variation. Fortunately for the range of 
VEV's considered here (i. e., large v 3 ), the only contributions to the parton level cross- 
sections turn out to be the diagrams which contain the Zj and H® propagators; the 
other terms are, in general, suppressed by several orders of magnitude.! Therefore 
the heavy lepton production cross-section is insensitive to variations in m pQ and 
m H ± . Here the P° mass will be set to 200GeV A . The corresponding H mass was 



chosen to be 215GeV for figure 4.3 and 212GeV for figures EO and |4.5|, which lies 



4 i.e., in the large 7J3 limit the couplings P°L + L~ — * 0, via equations C.18 and p. 38 , and 
H®L + L~ — * —(m L /v^) S^i , via equations |C. 17 and C.33, and equation 4.12 of Hewett and Rizzo jl3| 
for the C/3i's in this limit, to 0(l/i>3). 
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within the allowed regions on the m^o contour plots. Based on the very limited 
experimental constraints that do exist for supersymmetric models ]51| these appear 
to be very conservative choices. They also lead to fairly reasonable values for the 
masses. 

The next parameters that need to be fixed are the soft terms. Exactly how these 



terms should behave at low energy is not clear |[49|| . At the moment, their behaviour is 
very model dependent and unless supersymmetric particles are found this situation 
will most likely remain so. Here the soft terms will be treated parametrically as 
function of a single parameter m s . In particular, the soft terms will be assumed 
to be degenerate, Mf pa Af pa m s , with the exception of the A and A terms which 
were fixed by selecting the m p0 and m H ± masses. The selection of the soft terms 
in this way, including A and A for m p0 H ± £0(l)TeV, typify the generic outcome 
for the sfermion masses of most SUSY-breaking models, cf. JO], [HJ. In these models 
O{0.2)TeV < m s < O(10)TeV |^|. How low m s can be pushed depends upon the 
choice of VEV's (t>3 in particular, since it is relatively large). For the VEV's used here 
it was found m s > 0(400 — 450)GeV. In general the sfermions with light spartners 
have masses •> 0{m s ) , which are roughly degenerate (within O(b0)GeV) with their 
mass-eigenstate partners. The stops, t[ 2 , and the exotic squarks, q[ 2 , have splittings 
> 0(\ mt,q') , for fermion masses 0(200) <, m i)9 / < O(600)GeV, for low values of m s . 
As m s approaches 0(l)TeV all of the sfermion mass become degenerate and pa 
0(m s ) . 

Finally there is the matter of fixing the heavy fermion masses. The heavy quark 
masses, m q i , will be assumed degenerate as will the heavy charged lepton masses, 



mi' , such that m q i t i> > O[0A)TeV [51], p2| . The e' will be designated to play the 
role of the heavy charged leptons, L ± . 

Figures |4.6| through [4. 12| show the rapidity distribution at rapidity y=0 for pp — > 
gg — > L + L~ as a function of m L , for various scenarios. The rapidity, y, is defined 
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Figure 4.3: A plot of the mass contour lines as a function of m p „ and m H ±, for 
v i/ v 2 = 0.02 and v 3 /v 2 = 6.7 (m Z2 w A96GeV). The dashed curve in the upper left-hand 
corner is a plot of the zero of the Higgs potential above which it becomes positive. 
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Figure 4.4: A plot of the H® mass contour lines as a function of m p0 and m H ± , for V1/V2 = 
and v 3 /v 2 = 7.7 (m Z2 w 509Gey). The dashed curve in the upper left-hand corner is a 
of the zero of the Higgs potential above which it becomes positive. 
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Figure 4.5: A plot of the H® mass contour lines as a function of m p0 and m H ± , for V1/V2 = 
and v 3 /v 2 — 9.1 {m Z2 w A99GeV). The dashed curve in the upper left-hand corner is a 
of the zero of the Higgs potential above which it becomes positive. 
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Figure 4.6: Rapidity distribution at y = for heavy charged lepton production at the 
LHC (14TeV) as a function of heavy lepton mass, where vi/v 2 — 0.02, vy,/v 2 = 6.7, and 
m s — 400GeV. The mass spectrum for the non-SM particles involved in these processes 
are, m z ^ w A96GeV (T Z2 w 20.9GeV), m p0 w 200Gey (r p0 w 16.4GeV), m ff± w 215GeV, 
m H? « 94.3Gey (r ff? « 7.50 x 10- 3 GeV), m ff c « 200Gey (I^o w 16.5GeV), m ff c « 495Ge^ 
(T H0 w 0.230GeV), m g , = 200GeV. 
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Figure 4.7: Rapidity distribution at y = for heavy charged lepton production at the 
LHC (lATeV) as a function of heavy lepton mass, where v\/v2 — 0.5, v^/v2 = 7.7, and 
m s = AOOGeV. The mass spectrum for the non-SM particles involved in these processes are, 
m Z2 w 509GeV {Y z ^ w 21.5GeV), m p0 w 200GeV (r p0 w 2.52 x lO^GeV), m ff± w 212Gey, 
m H o « 75.4GeV (r ff? ~ 3.65 x 10- 3 GeV), m ff c « 212GeV (r ff0 w 7.49 x \Q- 2 GeV), 
m H0 w 507GeV (r ff0 « 0.198GeV), m ? , = 200GeV. 
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Figure 4.8: Rapidity distribution at y = for heavy charged lepton production at the 
LHC (lATeV) as a function of heavy lepton mass, where v\/v2 — 0.9, v^/v2 — 9.1, and 
m s = AOOGeV. The mass spectrum for the non-SM particles involved in these processes are, 
m Z2 w A99GeV {Y z ^ w 20.8GeV), m p0 w 200GeV (T p0 w 8.13 x 10- 3 GeV), m ff± w 212GeV, 
m H? « 52.3GeV (r ff? ~ 1-87 x lO~ 3 GeV), m H o ~ 216GeV (T H „ w 1.37 x l(T 2 GeV), 
m ff0 « 498GeV (T ffC w 0.130GeV), m ? , = 200GeV. 
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Figure 4.9: Rapidity distribution at y = for heavy charged lepton production at the 
LHC (14TeV) as a function of heavy lepton mass, where vi/v 2 — 0.02, v 3 /v 2 — 9.5, and 
m s — AbOGeV. The mass spectrum for the non-SM particles involved in these processes 
are, m Z2 w IQQGeV (T Z2 w 31.9GeV), m p0 w 200GeV {T p0 w 16.5GeV), m ff± w 215GeV, 
m H? « 94.6Gey (T H o « 7.49 x 10- 3 GeV), m ff c « 200Gey (T H o w 16.5GeV), m ff c « 700Gey 
(r ff0 w 1.04GeV), m ? , = 200G e y. 
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Figure 4.10: Rapidity distribution at y = for heavy charged lepton production at the 
LHC (14TeV) as a function of heavy lepton mass, where v\/v2 — 0.02, V3/U2 = 6.7, and 
m s — 400GeV. The mass spectrum is for the non-SM particles involved in these processes 
are, m Z2 w 496GeV (T Z2 w 19 AGeV), m p0 w 2QQGeV (T p0 w 16.4GeV), m ff± w 215GeV, 
m fl „ w M.SGeV (T H , w 7.50 x 10- 3 GeV), m H <, ~ 200Gey (T H o w 16.5GeV), m ff? « 495Ge^ 
(r ff0 w 0.138GeV), m 9 , = 600GeV. 
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Figure 4.11: Rapidity distribution at y = for heavy charged lepton production at the 
LHC (14TeV) as a function of heavy lepton mass, where v\/v2 = 0.02, v 3 /v 2 = 6.7, and 
m s — ITeV. The mass spectrum for the non-SM particles involved in these processes are, 
m Za w A96GeV (T Z2 w 20.9GeV), m p0 w 200Ge^ (r p0 w 16.4GeV), m ff± w 215GeV, 
m fl o w 94.3Gey (I^o w 7.50 x lO^GeV), m H <, ~ 200Gey (T H „ « IQ.hGeV), m H „ w 495Ge^ 
(r ff0 « 0.230GeV), m,, = 200GeV. 
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Figure 4.12: Rapidity distribution at y = for heavy charged lepton production at the 
Tevatron (1.8TeV) as a function of heavy lepton mass, where vi/v 2 = 0.02, v 3 /v 2 = 6.7, 
and m s = AQOGeV. The mass spectrum for the non-SM particles involved in these processes 
are, m Z2 w A96GeV (T Z2 w 20.9GeV), m p0 w 2QQGeV (T p0 w 16.4GeV), m ff± w 215GeV, 
m H o « 94.2Gey (r ff? « 7.50x 10- 3 GeV), m H o ~ 200Gey {T H ^ w 16.5GeV), m ff c « 495Ge^ 
(r ff0 w 0.230GeV), m ? , = 200GeV. 
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by 



y 



1, E + Pz 
- In ■ 
2 



(4.100) 



E- Pz 

and is related to the motion of the center of mass (cm.) with respect to the lab 
frame, where E is the cm. energy and p z is the net momentum of the incoming 
hadrons. The rapidity distribution, da/dy, is related to the parton level cross- 
section, a(gg — > L + L~), through 
da - M 



dy 



drG{^e\ Q 2 )G{^e~ y , Q 2 )a(T.s) 



(4.101) 



where r = s/s, r min = 4m 2 L /s, a/s is the center-of-mass energy, and G(x, Q 2 ) is 
the gluon structure function. The values of y/s have been set to 14TeV (LHC) 
for figures fO| through [4.1 1| , and l.STeV^ (Tevatron) for figure |4.12j . In these 
figures the SM couplings and masses where extracted from the PDG except 



for m t « 174 ± 13GeV [50]. For G(x,Q ) the leading order Duke and Owens 1.1 



(DO 1.1) |53], [54| gluon distribution was used. The results were compared with 



the next to leading order MRS A |55[ gluon distribution function, which yielded a 
negligible difference []52]. Although these results include squark mixing it was found 
that there was no significant change if mixing is not included |52[]. Since da/dy is flat 
about y = the relationship between ^\ y =o and the total cross-section is immediate. 
Therefore the total event rate for the pp — > gg — ► L + L~ production mechanism can 
be estimated from y = , i.e., 



a 



-i-dy w - ln(r„ 



(T 



(4.102) 



Figures |4.6| through |4.^ show d ' T 



dy \y=° 



for different VEV's ratios along the m Za 



O(500)Gel / contour line of figure [4.2| . Notice that as fi/f2 becomes comparable 
to V3/V2 the large V3 limit breaks down and the generally small qP° term starts to 
contribute (the q{Z\,2 — P°) contribution also grows quite significantly but remains 
a negligible contribution). Therefore for relatively large values of V1/V2 variations 
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in rripo (pa ?t%± up to at least 0(l)TeV) become important. However it is more 
natural to assume that the intra-generational Yukawa couplings are of the same order 
of magnitude and therefore for V1/V2 to be small. For the rest of these figures then, 
it will be assumed that Vi/v% = 0.02 . Figure O is the figure with the default values. 

Figure £0| shows what happens when a larger Z 2 mass of O(700GeV) (i.e., 
V3/V2 = 9.5) is used. For this figure m s had to be pushed up slightly to 450GeV ', in 
order to produce physical squark masses. The noticeable difference between this and 
all of the other figures is that the peak has broadened. This is expected since the 
Z2 can remain on-shell for larger values of m L . Notice that the H® resonance cut 
off seems to follow the Z 2 's. More precisely m H0 pa m Z2 for large V3. This becomes 
immediately evident when taking the large v 3 limits of the Hf and Zi mass mixing 
matrices, equations |4.30| and [4.47| respectively: i.e., 

JS8. - JS. '< - 1 «* * * f -i • < 4 - 103 ) 

which is in fairly good agreement with all of the figures. Also, the overall production 
is slightly suppressed due to the smallness of the gluon distribution function at large 
momentum fraction. 



Figure |4.10| shows what happens when the heavy quark mass was pushed up to 
QQOGeV . The effect is quite dramatic. To see why this is so, notice the slight kink 
in the curve around m L pa QOOGeV. There is also a much more significant kink in all 
of the other graphs around 200GeV, i.e., around m L ~ m ql . Further examination of 
the parton level cross-section shows that kink occurs when the heavy quarks in the 
loops can no longer be on shell. 



In figure |4.11| the scalar mass was pushed up to ITeV. Increasing m s has caused 
the terms involving the squarks to be supressed by several orders of magnitude. This 
supression occurs since there is not as much gluon luminosity available to allow these 
heavy squarks in the loops to be on shell, as there was for the light squarks. The qH® 
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term now enhances L + L production, below the m H0 threshold, as the destructive 
interference with qHf term, i.e., (q — q)H® , has been suppressed. 

Finally figure |4.12| shows what happens at y/s = l.&TeV, the Tevatron. The 



overall topology is the same as depicted in figure but the L + L~ production rate 
is dramatically reduced: very little gluon luminosity is available to produce these 
heavy particles. 

4.4 Discussion 



Figure [4. 13| gives a summary of the total cross-section for L + L production for the 



Eg model parameter space studied in the previous section. At the Tevatron it is 
immediately obvious that nothing will be seen, i.e., < 0(0.1) events/yr. For the 
LUC , over a reasonable range of parameter space, O(10 4±1 ) events/yr are expected. 
The MSSM results of Cieza Montalvo, et al, (CM) p| predicts O(10 5 ) events/yr, 
a factor of at least 10 more than our results |49], [52| . This is a rather surprising result 
since it was expected that the E$ event rate would be enhanced due to the greater 
number of heavy particles running around in the loop. The parameters in each 
model were varied to try and determine the difference. Unfortunately, it turns out 
that the E 6 parameters space suppresses L + L~ production, since v% is fairly large. 
This restriction causes the production to occur mainly through the Z 12 and H® 
terms. The MSSM has two neutral Higgses and one pseudo-scalar Higgs that are 
allowed to contribute to the processes. A very simple test of the Eg model was done 
by varying V3/V2, about t>i/t>2 = 0.2, that showed for V3/V2 = 2.8 and < lOOGeV^ 
a factor of 10 increase was obtained. However, this region of E 6 parameter space is 
forbidden, see figure pT2". 



The L ± — > u L W ± , v L H ± decay modes are expected to be similar for both mod- 
els, as these are SM-like decays. These modes depend upon the mass difference 
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Figure 4.13: Summary plot of results for the total L + L production cross-section at the 
LHC (y/s = UTeV, C - ltfpb^/yr) and the Tevatron (y/s = 1.8TeV, C - 10 2 pb- 1 /yr) 



energies as a function of m L . The hatched region is the results for figures 4.6 through 4.11 
and the dash line is the result for figure 4.12| . 
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A = m L — m u . For A < m w << m H ± the decays modes will be by virtual 
Ws, W* — > ff, and on shell for A > m w . Leptonic decays of the Ws offer 
the possibility of L^ 1 detection by measuring production with missing trans- 
verse momentum, j> T [16]. The SM backgrounds competing with these processes 
are pp — > t + t~ , W + W~ , Z°Z°. Studies have shown, using SM couplings, that it 
is possible to pull the L + L~ signals from the background for A > m w given suffi- 
ciently large event rates; it is much more difficult for A < m w [[U| |56], 0. Since 
the MSSM event rate is higher than in Eg , detection would more likely indicate 
MSSM candidates. If m H ± « 0{m w ) , then H ± — > fofj dominates, for naturally 
large values of tan (3. Since the Higgs likes to couple to massive particles, this would 
lead to multiple heavy jet events which in general would be very difficult to pull out 
of the background in either the MSSM or E 6 . Similar processes are expected to 
occur for the cases m w < A < m H ± and m w , m H± < A . For large enough m H± 
the sfermion channels also open, i.e., H ± —>■ fif* (e.g. ud*). The sfermions would 
eventually decay leaving only the lightest supersymmetric particles (LSP's), which 
will escape undetected along with the z/ L 's leaving lots of j> T . In fact, all of the 
aforementioned processes will lead to events with ^> T , as the z/ L 's will pass through 
the detector. 

In certain regions of the MSSM and E$ model parameter spaces it may be 
possible to distinguished between the two models. If L + L~ event rates are larger 
than those predicted by E 6 , then the likely candidate is the MSSM. Unlike the 
MSSM, it is possible for m H ± < m w in E$ [13] and therefore if H s are found 
in this mass range the more likely candidate would be Eg . Another possible way of 
telling the models apart is to look for sfermion production, i. e., L — > //* , // , which 
is unique to Eg , since L has opposite R parity to the other SM-like fermions in the 
27's (figure |1.5| ). The sfermion would eventually decay to an LSP which is stable 
(assuming R parity conservation), yielding jets+^> T , in general. Whether or not it is 
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possible to distinguish them from the MSSM and the SM-like backgrounds would 
require a much more detailed study, as the allowable parameter space for sfermion 
masses and Yukawa couplings is quite large. Finally, the MSSM does have fairly 
stringent unitarity constraints on the heavy lepton and heavy quark masses as a 
function of tan (3 fL6|, in particular m L < 1200 cos (3 . Therefore, it should be possible 
to eliminate (m L , tan f3) regions in the MSSM L~L + production cross-section plots, 
as a function of m L , such that only E 6 models are allowed. For example, assuming 
m H o > O{m0)GeV rules out the MSSM for m L >2A2GeV and tan/5> 5 [§2J. In 
the allowed MSSM region this gives an upper limit on the L + L~ production cross- 
section of O(10)pb p2[ , at the LHC . Also, in the MSSM there are phenomenological 
constraints on tan (3 which could allow for further restrictions. A more detailed 
study of these constraints has not been carried out. 

In closing, it should be pointed out that only a simple model of E 6 has been 
considered. It is possible for other E 6 models to produce results similar to the 
model studied here or to the MSSM . Therefore, in general, L + L~ production by 
gluon-gluon fusion should not be considered a definitive means of separating out the 
different models; several experiments would be required. 



4.5 Conclusions 

The pp — > gg — > L + L~ production cross-section was computed for a simple rank-5 
E 6 model. For a fairly conservative survey of the various parameters in the model we 
expect O(10 4±1 ) events /yr at the LHC; no events are expected at the Tevatron. 
It was found that the results were at least a factor of 10 less than the event rates 
predicted for the MSSM, due to the CDF and Dfi soft limits (i.e., assuming SM 



couplings) placed on m Z2 [50|. These constraints resulted in the H® 2 and P contri- 



butions to the L + L production rate being suppressed leaving only the H® and Z 1}2 
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to contribute. For certain regions in the MSSM and E 6 parameter spaces it was 
demonstrated that it is possible to distinguish between the two models, in principle. 
However, it should be pointed out that there are many candidate E 6 models which 
could yield overlapping results. 



Chapter 5 
Conclusions 



The chapter contains a summary discussion of chapters 2 through 4 followed by some 
closing remarks. 

5.1 Summary Discussions 

5.1.1 String-Flip and Flux-Bubble Models 

Chapters 2 and 3 of this thesis investigated various extensions and aspects of string- 
flip potential models in an attempt to extend these viable models of nuclear matter 
to obtain nuclear binding. 

In chapter 2 the SUg(3) and SUh{3) three-quark potential models were examined 
as possible alternatives to earlier SUt(2), SUh(2), and SUh(3) (these models included 
large multi-quark structures) models |4], |5| that may produce nuclear binding. Similar 
results to the SUe(2) models were obtained, but with no nuclear binding. It was 
found that generic features of all of the aforementioned models were: 

• nucleon gas at low densities with no van der Waals forces 
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• nucleon swelling and saturation of nuclear forces with increasing density 

• fermi gas of quarks at higher densities 

After a detailed analysis it was concluded that these models would have to be ex- 
tended to include perturbative QCD interactions between the quarks in order to 
achieve nuclear binding. Furthermore, it was hypothesized that the colour-hyperfine 
interactions would prove to be the most significant, as was indicated by some ear- 
lier work done with linked-cluster-expansion models [I7|. Chapter 3 marked the 
beginning of these extensions. 

In chapter 3, a new class of potential models, called the "flux-bubble models," 
were constructed that extended the old string-flip potential models to include per- 
turbative interactions. Some initial work was done using a flux-bubble model with 



SUg(2) by introducing colour- Coulomb interactions [58[|. Unfortunately, many prob- 
lems arose with trying to construct a new wave function which would take into 
account the local nature of the perturbative colour- Coulomb interactions. Also, the 
fact that we were dealing with a multi-quark system which required vast amounts 
of CPU resources only exacerbated the problem. As a result, a toy model of a 
mesonic-molecule (Q2) |0 was proposed, which consisted of two massive quarks 



(Q's) bound by an effective potential produced by two light quarks (g's), in order to 
sort out these difficulties. 

The Q2 system was analyzed using the old SUe(2) variational pseudo-hydrogen 
wave function || , \t> h , for the string-flip and flux-bubble models. The string-flip 
models were studied with and without a colour- Coulomb term added to the linear 
part of the potential, and the flux-bubble models were studied using SUi{2) (with a 
colour-Coulomb extension) with and without the colour-charge being allowed to flow 
from quark to quark. In general, it was found that the old \1/ h wave function pro- 
duced a very shallow Q2 potential well (0(3)MeV deep), not enough to bind the Q2 
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system. The exception was the SUe(2) model in which the colour-charge was allowed 
to flow. Here the system dissociated into two mesons, one heavy and one light. Next a 
new variational wave function was proposed, coined the "pseudo-hydrogen-molecular 
wave function," \& # 2 (i.e., its idea was motivated by the hydrogen-molecular wave 
function), in an attempt to take into account the local nature of the flux-bubble in- 
teractions. The aforementioned studies of the string-flip and flux-bubble models were 
repeated using ^ H , 2 . Surprisingly, the \l/# 2 wave function produced a well that was 
in general O(80)MeV deep, deep enough to bind quarks with masses m q > 0(m c ) . 
Here the exception to the case, SUi(2) with flowing colour- charge, did not dissociate. 
Although this produced a state higher in energy than with the *Sf h wave function, it 
demonstrated that the ^ h 2 wave function had the desired properties of localizing the 
colour-Coulomb interactions, as it bound the light quarks to the heavy quarks which 
would have otherwise floated away.0 Furthermore, the fact that the wells produced 
were, in general, much deeper seems to confirm earlier beliefs that the old SU((2) 
wave function used to model nuclear matter || was a poor choice |§ . 

The results of the analysis produced a more detailed understanding of intricacies 
involved in modeling nuclear matter. The situation as it now stands consists of 



several scenarios that were discussed in § 3.5 which may lead to nuclear binding. Of 



the various scenarios discussed it is believed that the most likely ones are, 

• the SUe(2) string-flip model using the many-body wave function, 

* ~ Permit ^Slater, 

• the SU C (3) flux-bubble model, with flowing colour-charge, using \1>, 



1 It should be re-emphasized, as stated earlier in chapter 3, that the SUg(2) with flowing colour- 
charge is not a physically viable model of nuclear matter. A more physical model would use 
SU C (3) and assume that heavy quarks are composites of two light quarks, in which case it becomes 
impossible for the system to dissociate. 
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both for a simple qq nucleoli model. These were the models that produced the most 
dramatic effects in the Q 2 system; the later model being only studied for SU C (2). 
Finally, it should be pointed out that this does not necessarily rule out the earlier 
hypotheses that the colour-hyperfme interactions are important, further investigation 
would be required. 



5.1.2 L + L Production 

In chapter 4, the pp — > gg — > L + L~ production cross-section was computed for 
a very simple (i.e., rank-5) superstring inspired Eq model. The L + L~ produc- 
tion rates were studied over the Eq model parameter space for, 0.02 < vi/v 2 < 
0.9, 6.7 < v 3 /v 2 < 9.1, 0.4 TeV < m s < ITeV , 200 GeV < m Q < 600 GeV , 
m p0 « 0(214) GeV , and m H± w O(200) GeV , as a function of m L for 0.05 TeV < 



m L < O(0.6 — 1) TeV . For this survey of the parameter space, we [^, |52| expect 
O(10 4±1 ) events/yr at the LHC and no events at the Tevatron. It was also found 
that the results were at least a factor of 10 less than the event rates predicted for 
the MSSM, due to the CDF and T>0 soft limits placed on m Z2 |jU| (i.e., v 3 /v 2 ), 
requiring m Z2 ^ O(500) GeV. These soft limits resulted in the if° 2 and P° contribu- 
tions to the total L + L~ production cross-section to become suppressed. It was also 
demonstrated, in principle, that for certain regions in the MSSM and E 6 parameter 
spaces it is possible to distinguish between the two models. Finally some discussion 
about the possible L decay modes ensued, with the general outcome resulting in 
jets + i> T . For these processes it is, in general, expected to be difficult to pull out 
signals from SM backgrounds; further study would be required. 
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5.2 Closing Remarks 

In closing, the string-flip potential models have not been ruled out as candidate 
models for nuclear matter. They correctly predict most of the bulk properties of nu- 
clear matter with the exception of nuclear binding. It is expected that the string-flip 
potential model with the inclusion of the newly proposed variational wave function, 
^ ~ Perml^^l $siater , along with perhaps some variant of the flux-bubble model, 
should produce nuclear binding. Finally, for L + L~ production (9(10 4±1 ) events/yr 
are expect at the LHC. If heavy leptons are discovered it may indicate the possi- 
ble existence of an E 6 model. However, in general, different experiments would be 
required to truly validate the model. 



Appendix A 

Three-quark K.E. Computations 



Given the wave function 

*(a,P,p) = X (<x,PMp), (A.l) 

where 

X(a,/3)=e ^ (A.2) 

and 

*<J>)= II l $ s c (p)l (A.3) 

ce{rgb} 

are the correlation and Fermi parts respectively, the kinetic energy can be split up 
as follows: into a correlation piece T c = 2T C _ S — Fq , a Fermi piece T F = 2T F _ s — F§ , 
and a mixed correlation-Fermi piece Tcf = — 2Fcf , where 

fc, - ^E«lox)), (A.4) 

^ - I^OT ■■•)>, (A.5) 

*2 = ^r<E(V,tax-) 2 ). (A.6) 

= ^T(E(V,ln$) 2 ), (A.7) 
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and 



^cf = <£(V 9 ln X ) • (V ? ln<t>)> . (A.8) 



The correlation pieces are straightforward to evaluate, and give the following 

{rgb} 

a 



- o f 3 " E Cgb 2 [(?r g - r hr )5 r£ + (r gb - f rg )5 g£ + (r br - r gb )<5 b£ ] 



^ {rgb} 



(A.9) 



with £ rgb = ^^r2 g + r 2 b + rg r , which yields 



and 



a 2 



XXV, In X) 2 = " ^ /? 2a E Cb" 4 [(% - ? br) 2 + (? gb - ? rg ) 2 + (r br - r gb ) 2 ] 



9 

t ^ {rgb} 



-«TE4" 2 . (A.10) 

{rgb} 



£ £ {rgb} 

•[(r rg — r br )<5rf + (r gb — r rg )<5 g ^ + (r br — r gb )5 b £] 
= - f F* E Cgb[(?r g - ? br ) 2 + (r gb - r rg ) 2 + (r br - r gb ) 2 ] 

{rgb} 

+18}Cgb 2 

= -a( a + 4)j3° E Cgb 2 , (A.11) 

{rgb} 



which imply 



and 



« {rgb} 



^ = ^-E«r 2 > (A.i3) 

MH 1 rgb 
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thus giving the desired result, i.e., 

T C = 2T C _ S = qr b 2 [«(l - (^ rgb ) Q ) + 4] . (A.14) 

Zm 1 {rgb} 

The Fermi pieces are also straight forward once the following identities are real- 
ized §: 

V £ ln|$ s J = y,o u V,o,j. (A.15) 

ran 

V, 2 ln|$ Sc | = E[^V^-E^-(V £ 0L)-(V^)0n,], (A.17) 

ij ran 

where 0^ = (0 T )« , and 0^ = sin(^ h*j • ?j + <5j) . These imply 

27T 

Vi&j = —niSji^ij (A. 18) 



and 



47T 2 

V^ = -^4, (A. 19) 



where 0^- = (j>iifj + j^-Si) • Therefore 

47T 2 

E Vf ln I $ s c I = - -jj E 0« Wi + E M'A ■ "jo'y.oy,. (A.20) 

£ ij ijk 



and 



imply 



4tt 2 

E(V,ln |$ Sc |) 2 = — E Mtfk ■ n.,0',,0,, , (A.21) 



T (c) _ 9T (c) p2 ( c ) _ 2?r I =*2iT 



Tr(0 T iV 2 0) = -^-Tr(Af c 2 ), (A.22) 



m q L 2 VY " ' m q L 2 
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where (N^)ij = n?6y . Thus the desired result is obtained by summing over the 
quark colour degrees of freedom, i.e., J2 c Tf \ which implies 

Tr = $? 3 '- (A - 23) 

Finally, the mixed correlation- Fermi result, given by equation ( |2.25| ), is simply 
obtained, via equation ( |A.8|) , by taking the inner product of equations (|A.9| ) and 
(|A.15|) , summing over i and q, and using the identity given by equation (|A.18| ). 



Appendix B 
Algorithms 



This appendix contains useful algorithms that were used in some of the Monte Carlo 
routines in this thesis. 

B.l Generating Slater Wave Functions 

This section shows how to generate the Slater wave function for a given number of 
fermions, rip. 

The first thing to notice is that equation ( |2.13| ) can be rewritten as, 




where 



77(h) 



3e(ni) + 2e(n 2 ) + e(n 3 ) , 



(B.2) 



is the signature function, and 



-1 



if x < 



e(x) = { 



if x = 



(B.3) 



+ 1 



if x > 
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This definition is much more useful for programming purposes because it gives a 
unique wave function for a given lattice vector, n. 

The final step is to develop an algorithm to generate the lattice vectors: consider 
a set of lattice vectors, n, whose components are given by 

n a = -k,...,k, (B.4) 
where k is an integer. Then the total number of states is (2k + l) 3 , or 

fc = -r# 3 -l. (B.5) 
Therefore the algorithm is as follows: 

• Generate (2k + l) 3 vectors n, with components n a = —k, . . . , k. 

• Sort the vectors in increasing order of magnitude squared, n 2 . 

The sorting step makes sure that the states, 0g(r) , are increasing in energy: since 

EFermi ~ |3| 2 . 



B.2 Distributed Minimization 

The section gives a detailed outline of a distributed minimization algorithm that was 
used for the work done in chapter 

Consider a <i-dimensional hypersurface that is being fitted locally to a paraboloid 



as depicted in figure |3.4|, then the minimum number of points required for the fitting 



is 



p= -d(d - 1) + 2d + 1 , (B.6) 



and the distributed minimization algorithm, outlined in § |3.4.1|, is as follows: 



1. Submit k > 0(2p) + ip "sample" points, x K , s.t. ip ~ 0(m) . 
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2. Wait for f(x v ) points to arrive, s.t. rj ~ k — ip . 

3. Push the results onto stack S. 

4. Fit each point f(xk), k = 1, ... ,r], using their N > 0(2p) nearest neighbor 
points in S, where N is chosen s.t. xl is minimized. 

5. Fit S in its entirety. 

6. Push the critical points, x p , corresponding to parabolic minima onto stack M. 

7. Purge any redundancy out of M. 

8. Submit the newly found critical points x p eM that survived the purging. 

9. If p = and M is not empty choose a point in M with the lowest xl an d refit 
it to obtain a new critical point. 

10. If this point is a minimum and is not in M then push a copy of it onto M and 
submit it as x p= i . 

11. If this is the first pass submit (ip — p)6(ip — p) "sample" points, otherwise if 
p = and there are tp — 1 points pending submit a "sample" point. 

12. Repeat all of the above steps, except steps 1 & 5, with rj — 1 until some 
convergence criterion is meet. 

13. Flush all the computers. 

14. Return the lowest result, (f(x m i n ),x min ) , off of stack S. 

Step 1 is the initialization step which consists of overloading the system, of m com- 
puters. In step 2 of the initialization phase the routine waits for 0(2p) points to 
arrive. Once these points have arrived the system is left occupied processing the 
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remaining ip points. Therefore ip is chosen by the user, who decides how many of 
the m computers they want to use (keep loaded). Although all that is required for 
a fitting to a <i-dimensional paraboloid is p points it was found that 0(2p) points 
gave much better statistics. In steps 2 and 3 all of the incoming Monte Carlo are 
fitted to a d-dimensional paraboloid and stored in a data base: i.e., stack S. As this 
routine passes through several iterations the data base, S, builds and the statistics 
becomes better. To take advantage of this as many points as possible are fitted in 
the neighborhood of a newly arrived data point, Xk ■ The optimum number of points 
is determined by minimizing xl such that there are at least 0(2p) of them. Step 5 
is an optional step which attempts to fit the entire contents of S. This step is only 
done in the initialization phase of the routine. Steps 6 through 8 submit only the 
critical points that correspond to a parabolic minimum and that have not already 
been submitted. Therefore stack M contains a list of points (along with other per- 
tinent fitting information) that have already been submitted. If these steps fail to 
submit any points, steps 9 through 11 take over. Steps 9 and 10 attempt to find a 
new parabolic minima to submit by refitting about the neighborhood of a point in 
M, using the data in S, that has the lowest xl ■ If these steps fail then step 11 takes 
over by submitting more "sample" points, as follows. If on the first pass no points 
were submitted the system becomes doubly loaded. On future passes no points are 
submitted until the overloading has vanished, at which point the system is kept 
optimally loaded. Step 12 checks its convergence criterion to determine whether or 
not the routine should continue. If the convergence criterion is met, steps 13 and 14 
clean up the system and returns the result, which is taken as the lowest value in S. 
The "sample" points were generated by using a Sobol quasi-random sequence |f)9f 



that ensured a uniform sampling of the region being searched, such that each point 
was only visited once. Therefore if all of the attempts failed at fitting parabaloids, 
to predict the next step, the searching region would be continually sampled: i.e., if 
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becomes effectively a grid search. 

In Step 11 an option was used in which, under the condition that after the first 
pass if M was not empty the sampling was done locally about the point that failed 
to yield a minimum in step 6. The locality of the sampling was determined by 
the uncertainty in the position of critical point from its past parabolic fit, in the 
previous pass. Therefore, assuming several passes lead to execution of this step, 
more statistics would be gathered locally about this region and the search would 
either converge or wonder away. In order to make sure that this searching thread 
did not leave the box, certain boundary conditions were placed on it: i.e., if the 
point was near the edge of the box only intersection of its error box with the search 
region was sampled. 

For the work presented in this thesis, a convergence criterion was used in which 
the search would end if the <, 0(1) , for the point with the lowest xl i n M. Later 
work showed that minimization using a weighted variance produced more satisfactory 
results. This method reduced the computation time by ~ 1/2 . 



Appendix C 

Couplings and Widths for 

a(gg^L+L-) 

This appendix gives a summary of the calculations that were used to obtain the 



couplings and the widths for the cr(gg — ► L + L ) matrix elements given in § 4.3 



C.l The Couplings 

In this section the calculations of the vertex factors used to obtain the gg —>■ L + L~ 



matrix elements, given in § 4.3, are summarized. 



For the Z 12 exchange diagrams shown in figure |4.1| .a the following vertex factors 
were used 



, f 



-g 



yl — X\y 

f 



7 



»[C f L Zi P L + C f R Zi P R ), (C.l) 
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where i = 1, 2, 



= -(1 -75) 
^ = ^(1+75) 



(C.2) 
(C.3) 



and C q L Zi and Cj^ were the couplings used in equation 4.82 . The gauge- fermion 
interaction Lagrangian for SU(2)l <8> U(l)y <£> U(1)e is given by 



(C.4) 



where the ip^s are two-component spinors, see equation B.2 of Haber and Kane 
(HK) g§. Defining g = g L , g' = g Y , g" = g E , and F = Y Yx {= %Q - f 3 ), and using 
the identities 



aL^ + bY^ = (acos9w — bsm9w)Z fJi + (asinOw + bcos9w)A fl , (C.5) 



(C.6) 



then equation |C.4| becomes 
5 



1 



f ()s (T 3 - Q x w )Z, + - g"Y E Z'^ \ [^ fL)i a^ iMj + V(/ E)i ^(/ £ )J 



where T 3 = f 3 /2 , r» = , tan 6^ = g'/g , and rr^ = sin 2 . Noting that 



(T 3 -Q^M> 



-(T 3 - Qa^w)l/fl > > 
-F £ |/*>, 



(C.7) 

(C.8) 
(C.9) 



yields 

Ant 3 — 



— ^-(r 3 - Q + ~ g"Y E z' \ [i'(f L ) t (T fl ^ i f L)j - ij( fR )^(f R)j } . 

COS (7py Z J jj 

(CIO) 
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Using the following identities 

to convert from two- component to four-component spinor notation yields 



(C.11) 
(C.12) 



-9 



J2 f { ( T 3a ~ e f x w) % 



i (g' 1 



see equations A. 28 of HK Then the Z-Z'- vertex factor is 



(C.13) 




-g 



yl — X\y 



7 



nc f L z ' fz 'p L +c f R z ' fz 'p R } 



(C.14) 



where the C[ Z J Z 's are defined by equations [4.84| and [4.85f Using the inverse of 



transformations fO and 4.2 



Z' 1 / oos r; — sun 
Z / \ sm0 cos( 



^ ^2 



(C.15) 



yields the desired result 



-g 



Y. £ JZiC f A Zi P A f 



(C.16) 



v 7 ! - j =1 A=L>R 

i.e., vertex factor CO . 

For the Hf 2 s an d -P exchange diagrams shown in figure |4.1| .b the following 
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2m w 



(C.17) 



and 




(C.18) 



respectively, where z = 1,2,3. The K^ H i and K^ p ° couplings are obtained from 
the Yukawa interaction part of the Lagrangian given by equation |4j]: noting that 
Sij = (ir 3 )ij and plugging W, equation £D| into equation [17| yields 



C 



Yuk 



D 



+ M[4>i{^d L ^di + tpditpd L ) + <t>i$d3d% + ^d-J)d L )\ 

+ A 5 [0° ^d> + i>d> i>d' c ) + 03* fe M + ipd' */V c )] } > 



(C.19) 



(C.20) 
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and similarly for the other generations. Defining 

f. 

/ = 



and using the following identities 



see equations A. 24, A. 25, and A. 28 of HK f45| , implies 



A-Pl/2 , 
/2-Pr/i , 



-\K [#/(l-7fc)/ + # /(I + 7s)/] 
-A, [Re(0°)/7-iIm(0;)/ T5 /], 

-73 A* (oy/-'4/",[). 



(C.21) 



(C.22) 
(C.23) 

(C.24) 
(C.25) 
(C.26) 
(C.27) 



Expanding the 0°'s in terms of their physical fields, equations |4.39| through |4.44j , 
yields 



Equation 4.54 

1 3 

— y=J2 {AiC/ 2 ., + C/y(A 2 + A 3 ee) + U 3j (X 4 ee' + A 5 d'tf)} #° 



IK 



+-^={Aii;i3 M7 5 n + w 23 (A 2 dj 5 d + A 3 e7 5 e) + v i 2 (A 4 e'7 5 e' + A 5 d'jsd'tyP , 

(C.28) 



where k = \j Jv\v\ + t> 2 t> 3 . The couplings can now be read directly and give, wo 



equations 4.55 through 4.59 



sm p 



(C.29) 
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K 



cos (3 



(C.30) 



K 



d'H° 



2m w 



u 3i , 



(C.31) 



K 



-.H° 



cos (3 



(C.32) 



7T 



for the scalar Higgs fields, if? , and 



5^3 



17: 



3i 



(C.33) 



A" 



sin {3 



KV 13 



(C.34) 



7\ 



cos (3 



(C.35) 



d'p° 



2m 



w 



gv3 



KV12 



(C.36) 



K 



cos /3 



(C.37) 



iT 



2m 



S^3 



(C.38) 



for pseudo-scalar Higgs fields, P° . 

For the H l23 exchange diagrams shown in figures |4.1| .c and |4j].d the following 
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vertex factors were used 



" ■ JS = --f^K^, (C.39) 

V 

off 

where A, 5 = L, R . The couplings, which were obtained directly by using 



Mathematica p6| to generate equation |4.8| and to extract the couplings from it, are 
as follows: 

4? = (y 2 - y 2 + y 2 ) vuvi + y 2 + \g 2 - y 2 - a?) u 2i v 2 

-t 8 g" 2 u 3i v 3 , (C.40) 



d'H? 



-y'U 3i u 3 , (C.41) 

K, 



\ [{U 3i pi + U u u 3 )\ - ^[2U 2l A u \\ x , (C.42) 



V LR — 2 

4? = (y 2 + \g 2 + y 2 -^)Uuvi + (lg' ,2 -\g 2 -y 2 )u 2l v 2 

-^g" 2 u 3l v 3) (C.43) 



V RR — K 36 



(y - y + \d u u v x - (§</' + §</ ) u 2l v 2 



+y 2 U 3i v 3 , (C.44) 
4? = \[{U 3l v 2 + U 2l v 3 )\- V2U u A d }\ 2 , (C.45) 

= -y 2 + y 2 )u li u 1 -(ig" 2 -y 2 )u 2i v 2 



+ (lg" - \f) U 3i v 3 , (C.46) 

4^ = -(^ ,/2 -|/)^i^i-(^ //2 + |/)f/ 2 ^ 2 

+ (|/ 2 -AD^3^3, (C47) 



4? = H( f/ 2^i + ^^)A- v^C/a^'jAs, (C.48) 
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for the squark-Higgs couplings, and 



eH° 
K LL 



"RR 



+ 3 \g" 2 u 3i u 3 , 



' 2 „// 



\g 2 -\g')u 2i u 2 



(±g" + \g' - A3) u u v, + (§</' 2 - §/) r/ a ^ 2 



K 



eH° 
LR 



\ [{U 3i v 2 + U 2i v 3 )\ - V2U u A e }\ 



2 ■, 



UeH° 

"RR 

VeHl 
V LR 



-ihg" + W + \g'l u u ^ - {\g ,r - \g 2 - \g r ) u x ^ 



+■^9" U 3i u 3 , 

^g"\Uu v\ + 4U X v 2 - 5U 3i 1/3) , 

0. 



(C.49) 

(C.50) 
(C.51) 

(C.52) 
(C.53) 
(C.54) 



~ fH° 

for the slepton-Higgs couplings. The mass eigenstate couplings K 1 2 1 , given by 
equations |4.89| and |4.90| , were obtained by substituting 

h I / cos6>? -sin Of \ fi 



f. 



R 



sin 1 



cos 6 




(C.55) 



which is just the inverse of equation |4.77| , into the scalar potential and extracting 
them via Mathematica |36 . 



The corresponding pseudo-scalar-Higgses couplings 



I K 



d'P° 
AB 



—I 



gm z 



\/l - x w 



K 



gP° 
AB 



(C.56) 



"fA 

are obtained in a similar fashion as above: i.e. 



h AB 



— G 



>\ + z/f)A + V2 A u v 3 cot p] Ai , (C.57) 
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Jp° 

V AB 



K 



d'P° 
AB 



~AB 



'AB 



"3 

2 A 2 2 + ^3 



; )A + A^ 3 tan/3 



A 2 



^3 



\5 • 



for squark-pseudo-Higgs couplings, and 



AB 



2 Jz/^ + z/V 3 2 



+ f|)A + V2 A e u 3 tanp 



A3 



K. 



£> e P° 

AB 







for the slepton-pseudo-Higgs couplings, where 

1 if A = L, B = R 



~AB 



if A = B 
-1 if A = R, B = L 



(C.58) 
(C.59) 

(C.60) 
(C.61) 



(C.62) 



fpO 

In general, these couplings will also mix to give the mass eigenstate couplings K{ 2 ; 



which are defined in a similar fashion to equations 4.8S and 4.90 



C.2 The Widths 



In this section all of the tree level two-body decay widths for Z 2 , Hf , and Pq are 
computed. The generic two body decay formula is given by |K) 



S\M ab \ 2 



lQ-Km 



fiab , 







m a 
O m 
m b 



(C.63) 



where rrii , i = 0, a, b , are the masses of the particles, pi , in the decay process 
Po -> PaPb , 

9 I'm 2 _|_ <m?l f r>)2 — m?"l2 

(C.64) 



/3a 



\ 



_ 2(mg + m 2 ) (mg - m 2 ) 2 
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such that p a b = fi a H a = b , S is & symmetry factor for the out going particles, p a 
and pb , and M. ab is the amplitude for the process. 



C.2.1 Yz 2 

For the Z 2 width the following processes need to be computed: 
Z 2 — > W+W-, Z,Hl W ± H^,q i q i , hk, x°iX% xfxj , QiQp U], H?H°, H+H~, P°H° . 



The Zo — > W + W width, which can be found in Hewett and Rizzo ITS 



is given 



by 



q 2 m 7n sin 2 /'m 7o \ , 

r z 2 -» w+w- = f—f \ ) p 3 w 

1927r(l — xw) V m z J 
The Z 2 — > /j/j vertex factors are given by 



1 + 20 



rn 



IV 



+ 12 



(C.65) 




(C.66) 



u f (p) 



which were obtained from equation |C.16| by converting to the V — A basis: i.e., 



aP L + bP R = v f - a/75 



(C.67) 



where 



Vf = \{a + b) 
a f = \{a-b) 



(C.68) 
(C.69) 



APPENDIX C. COUPLINGS AND WIDTHS FOR a{GG -> L + L~ 



146 



which yields 



v f 
a .f 



2yl — 
1 

2yl — 



(cf 2 + cf 2 ; 



(cf 2 -cr)- 



Therefore 



1 „2 
3 



/EE W)7>/ - a/75)«(p)e M (?, A)p 

A spin 

lg 2 [v}{m 2 Z2 + 2m})+a}{mj 



4 - 4mJ) 



Plugging this into equation |C.63| gives 

^2 



T(Z 2 -> ff)= Cf ^-m z J f 

where c/ is a colour factor which is 3 for quarks and 1 for leptons. 
The Z 2 — > qiq*, hi* vertex factors are given by 



m 



z 2 . 



Amy 



ni 



z 2 



(C.70) 
(C.71) 



(C.72) 



(C.73) 



/ ft 



Pi 



-igKijirfj - Pi) 



Pi 



where fk=i,2 are sfermion mass eigenstates, and 

Vf + cif cos 29 f ; if i = j = 1 



(C.74) 



(C.75) 



Vf — df cos 26*^ ; if z = j = 2 • 
— cif sin26*j ; if i j 
The vertex factor is obtained by followings steps similar to equations C.4 through 



D 



n — -E E h z? 

y/L-Xw i= XA=L,R 



(C.76) 
(C.77) 
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followed by transforming the sfermions to their mass eigenstates by using equa- 
tion p.55| , 



£ mt D 7 =JE= J2{[C f L Zi cos 2 6 f - + C f R Zt sin 2 Ofi ft gji 



yjl — X W i=1 



+ [Cf* sin 2 0,- + cos 2 9 f ] ft d,f 2 
~ § [C f L Zi - <5£*] sin 26 f (ft dj 2 + ft djx)} Z? , (C.78) 

and then changing to the V — A basis, via equations |U.67| through |C71| , to obtain 



C- mt D -ig £ KijftdJjZj;. (C.79) 

i,j,k=l 



Therefore 



V<«M,. (c.80) 



— 3i/ ""Za^ij^fif. 

Plugging this into equation |C.63| gives 



r(^-/^) = C/ 4^4%- (C8i) 

For the range of VEV's that will be considered here (i.e., large V3 in particular) 
the Z 2 -> Z r Hf, W ± H^ 1 H°H°, H + H~, P°H? widths can be approximated by 

where the H^H® are kinematically forbidden or suppressed |fl3f 



The Z 2 — > XiXj widths are quite complicated to compute, due to the 

complex nature of the mass matrices, and can contribute as much as 10-20% to the 
total width, neglecting phase space suppression [13|]. Here its contribution will be 
taken as 15%; this proved to have no noticeable impact on L + L~ production. 
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C.2.2 r^o 

For the iJf widths the following processes need to be computed: 



The H® — > ZjZk, W + W vertex factors are given by 




where 



h9 

^VaVb 9fMV 



(C.83) 



c 



Z\Z2 
Z2Z2 



Czz cos 2 — 2 Cz ' Z i sin 20 + sin 2 . 
C%* z , cos 20 - (Czz - Cz'z') sin 2 > 
C zz sin + 2 C z ^/ sin 20 + C z / Z / cos « 



(C.84) 
(C.85) 
(C.86) 



for the Zi's, with 



and 



C 



% 

z z 



^ZZ' 



^Z'Z' 



- (g cos 9 W + g sin 6 W ) 2 (U u vi + U 2i v 2 ) 



g" 



(g cos 9 W + g' sin 6 W ) (U u vi - &J 2i v 2 ) 



%- {Uuvi + lQU 2l v 2 + 25U 3l v 3 ) 
6b 



c w+w- = 7 (^1^1 + 



(C.87) 
(C.88) 

(C.89) 
(C.90) 



for the If' s. The vertex factors, Cy o V , were obtained, via Mathematica [ 36jl , by 
plugging equation |C15| and equations |4.39| through |4.44| into the kinetic terms for 
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the scalar-Higgs fields, equation |4.22|. Therefore 



E |e M (p a ,A )<V 6 <(p 6 ,A 6 )^| 2 



H° 



171 H? C VaV h * 



A(m a m b ) 2 
which yields, via equation |(J.6ci 



i _ 2(m 2 a + m 2 b ) (ml + m 2 b ) 2 + 8(m a m h ) 2 



(C.91) 



647r(m a mfc) 2 



_ 2(m 2 + m 2 ) (m 2 + m 2 b ) 2 + 8(m a m b f 

2 ' _,4 



(C.92) 



where 5 = | for identical Zj's, otherwise 5 = 1. 
The — > Zj/j decay width is 



r(i*<W/) = ^(^Yff /fl? V 

1 1 ^ 32vr ^m^y 



H9 



(C.93) 



via equation |(J.63| with amplitude 



m^p — 4m 2 ] , 



(C.94) 



where the K^ H i couplings are defined by equation p. 17 . 
For the scalar processes H® — > 00* the vertex factor is 



Pb/ / 

C,X , (C95) 



Pa \ \ ^ 



which yields the following decay width 
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via equation |C.63|, with amplitude 



(C.97) 



For Hf — > g^g^, Z^. the vertex factors are 



gm z K f~ H ° 



(C.98) 



where the K^ k i couplings are given by equations |C.29| through |C.33| . For H\ 



H?H%, H + H~, P°P° the vertex factors are: 



- XA (U 12 U 21 U 31 + U 1X U 22 U 31 + U n U 21 U 32 ) 



12 



^(g" 2 + 9g 2 + 9g' 2 )U 2 n -(±g" 



4 



+^ ,,2 ^i-^ 2 (c/| 1+ c/3 2 a ; 



E7- 



11 



-1 



-8 5 



1 



^^ 2 c/3it/ 32 - A 2 (ry 21 c/ 22 + C/31C/32) 

3o 



+ *7: 



'22 



^ <T + ^ + -/ ) Ul x - -(16<T + 9g 2 + 9/)^ 



1 

2' 



+ tf 



21 



+ ^g" 2 u 31 u 32 - x 2 (u n u 12 + u 31 u 32 ) 



v 2 



+ U; 



'31 



+ u 32 

1 

2 



3o 

-g" 2 (U 2 u + 4t£ - 15L&) - l -{U 2 n + fft)] } v 3 , 



72" v 11 Zi Ji/ 2 

A/4 (U 13 U 21 U 31 + UnUMx + U n U 21 U 33 ) 

1 / //2 . „ 9 . „ /2 . T T o / 1 „2 1 



(C.99) 



'13 



-far + ^ 2 + 9/)^ - - ^ 2 - j<r a 



24 



31/ 



17 



11 



+ 720" Eft - 5 A M + U > 

+ 4^ 2f/ 31^33 " A 2 (t/ 21 f/ 23 + U 31 U 33 ) 
60 



21 



-g-P + ^0 + 4^ ) ^21^3 



£7 



23 



^g" + ^ 2 + \g" 2 ) U 2 U - ^(16/ 2 + 9g 2 + 9g' 2 )U 2 21 
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H"0 



U H°H° 



U- 



21 



+ \g" 2 U 31 U 33 - X 2 (U n Un + U 31 U 33 )^ } t; 2 



'31 



g" 2 (U u U 13 + 4C/ 21 C/ 23 ) - A 2 (C/ n C/ 13 + C/21C/23) 



+ C/33 
1 



36 

^/ 2 (t/ 1 2 1 + 4C/ 2 2 1 -15C/Ii)-^(C/i 2 1 + C/ ; 



^3 



(C.100) 



\A (U 13 U 22 U 31 + U 12 U 23 U 31 + U 13 U 21 U 32 



+u u u 23 u 32 + u 12 u 21 u 33 + u n u 22 u 33 ) 

- { i^g" 2 - \g 2 - \g" 2 ) (u 12 u 21 + u u u 22 )u 23 

+U 13 [^(g" 2 + V + 9g' 2 )U u U 12 + - ^ 2 - i/) U 21 U, 



^ </ /2 £/ 3 iC/ 32 + A 2 (f/ 21 f/ 22 + U 31 U 32 ) 



- ^g" 2 iu 12 u 31 + u u u 32 )u 33 } Vl 

-\ 2 (U l2 U 21 U 23 + U n U 22 U 23 + U 12 U 31 U 33 + U n U 32 U 33 ) Vl 



1 „// 2 



4 



+£/ 



2:-! 



1 „2 



g<T - ^ - ^" 2 ) C/llC/l2 + (^" 2 + \ 9 2 + ^ ) I W 22 



\ g" 2 U 3l U 32 + \\U u U l2 + C/31C/32) 



- g tf" (C/22C/31 + U 21 U 32 )U 33 l V 2 

-A 2 (f/i 2 f/i 3 f/ 2 i + U U U 13 U 22 + U 22 U 31 U 33 + U 21 U 32 U 33 )v 2 



+ 
+ 
+ 



^ g"\u 12 U 13 + 4C/ 22 C/ 23 ) - A 2 (f/ 12 C/ 13 + C/ 22 C/ 23 ) 



■ g"\UnU 13 + AU 21 U 23 ) - X 2 (U n U 13 + C/ 2 iC/ 23 ) 



5 „j 
36 

|j / 2 (f/nf/i2 + 4C/ 21 C/ 22 - 15C/ 3 iC/ 32 ) 



C/31 
C/32 



-A 2 (t/iit/i2 + t/ 2 it/ 



22/ 



^33} «3, 



(C.101) 



A/4 (U 13 U 22 U 32 + U 12 U 23 U 32 + U l2 U 22 U 33 ) 
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+ U 



'13 



2 

22 



J. 2 



+ ^g" 2 u 32 u 33 - A 2 (f/ 22 f/ 23 + tfotfo) 



g-9" + 4/ + ) ^22^23 



■'23 



18 P 8^ 8^ J 12 _ 24^ P 9 9 ' 22 



+ ^X-^ 2 2 + [/ 3 2 2 ; 



22 



33 , 



+ 7:9" 2 U 32 U 33 - X 2 (U 12 U V3 + L/^ 
y 

+ [u 32 [^g"\u 12 u l3 + 4f/ 22 f/ 23 ) - X 2 (U 12 U 13 + c/ 22 c/ : 



-g-^ + ^g +^g ) u u u 13 



v 2 



33 



^"> 2 2 + 4£/ 2 2 - 15U 2 2 ) - i(t/ 2 2 + C/ 2 2 ) 



23 , 



V3 



(C.102) 



for the neutral-scalar-Higgses; 

,ffp -1 



C 



H+H- 



{^(3^ + / -AX 2 ){U u + U 2l cot (3) 



4(1 + cot 2 0) 
+ (V - g' 2 + ^ //2 ) (vi^k cot 2 + v 2 U 2i ) 
+y ,2 [vi{v 1 U li + 16U 2i cot/3) - 5v 3 (l + 4cot 2 /3)£/ 3l ] 
+4 (A 2 t> 3 (1 + cot 2 (3) + X cot [3)U 3i , 



(C.103) 



for the charged-scalar-Higgses; 



C 



H9 

pQpO 



2{ylvl + v 2 vl) 



2 0R 1 (g"\ 

m 7 cos 2p — m 

z nW 



2 \2 V 1 V 2 



w 



4 



(viU u + v 2 U 2i ) 



v x v 2i 



-X(vfU u + vp 2i ) - X 2 v 3 (vf + v 2 )U 3i - A/4 (vtUu + v 2 U 2i + v 3 U 3i ) 

v 3 



2„,2 



{ Vl U u + Av 2 U 2l - 5v 3 U 3i ) + v 3 {Avl + v 2 )f/ 3i 



, (C.104) 



for the pseudo-scalar-Higgses, which were all extracted, via Mathematica [BO], by 



plugging equations [4.39| through 4.44 for the physical Higgs fields into the Higgs 



potential, equation 4.18 
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The H® — > Xj'Xfc? XjXk decay processes are quite complicated to compute. Here 
a simple approximation was made in which for m p0 < O(500)GeV its contribution 



to the width was 15%, otherwise 50% [H, p2[ . This addition had a negligible affect 
on L + L~ production, since m p0 = 200GeV. 



C.2.3 Tpo 

For the P° width the following processes need to be computed: 



pi. 



ZiZj, W ± H^, qi q h hi, XiX°j, XiXj, QiTj, U] ■ 



For P° — > ZiZj, W ± H T the widths are zero since here m p0 < m Z2 and m p0 rs 
m H± : see figures [12] through [15] and discussion therein. 
The P° — > Zj/j decay widths are 

T(P° - ff) = %f (^X Kf p ° 2 f3 p0 m p0 , (C.105) 
32tt V m w// 

via equation |C63| , with amplitudes 

*W 

where the K^ p ° couplings defined by equation |C 18 
The P° — > (jjql, ljl k decay widths are 



m 2 P o , 



Cf9 



2 m| 



16tt(1 — x w )m p0 



via equation |C96|, with vertex factors 



Cf, 



9 ™> z R fpo 



where the Kj P couplings are given by equations |(J.34| through |(J.38 



(C.106) 



(C.107) 



(C.108) 



In this work m p0 was fixed at 200GeK. At this mass P° — > XiXji XiXj decays 
are suppressed 
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